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ABSTRACT 

Results  of  a  series  of  experiments  on  optically  illuminated 
superconducting  tunnel  Junctions  are  presented.  In  the  first  chapters, 
the  question  of  what  nonequllibrium  quasiparticle  energy  distribution 
results  from  optical  Illumination  on  a  thin  film  is  addressed.  A  new 
technique  using  tunneling  current-voltage  characteristics  is  developed 
which  accurately  measures  the  nonequllibrium  quasiparticle  occupation 
function,  f(E).  The  extremely  good  resolution  of  this  technique  is 
sufficient  to  compare  the  experimental  situation  to  contrasting 
theoretical  models  of  raised  temperature  (T*)  and  shifted 
electro-chemical  potential  Data  taken  on  the  Al-PbBl  tunnel 
junctions  illuminated  for  this  study  are  well  described  by  the  T*  model. 

Stiidy  of  Illuminated  tunnel  junctions  lead  to  the  discovery  of  a  new 
thermoelectric  effect  for  tunneling  through  oxide  barriers  between  metals 
at  different  temperatures.  The  final  chapters  are  devoted  to  examination 
of  this  thermoelectric  effect.  Two  types  of  experiments  were  performed 
to  test  for  the  existence  of  this  effect,  with  positive  results.  A 
self-consistent  Rothwarf -Taylor  model  is  presented  which  correctly 
predicts  the  sign,  magnitude,  temperature  dependence  and  laser-power 
dependence  of  the  thermoelectric  current.  The  thermoelectric  effect 
described  here  appears  to  resolve  long-standing  discrepancies  between 
experimental  resiilts  and  theoretical  predictions  for  a  series  of  point 
contact  experiments. 
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FIGURE  l.l.  I(V)  for  an  Al-PbBi  tunnel  junction  at  2.08  K.  Dashed 
lines  indicate  hysteretic  voltage  jumps.  The  current  is  roughly 
exponential  at  small  voltages. 

FIGURE  2.1.  Tunnel  junction  layout. 

FIGURE  2.2.  Substrate  holder  showing  electrical  contacts  in  upper 
block. 

FIGURE  2.3.  Block  dlaigram  of  LSI-11/2  microcomputer  system. 

FIGURE  2.4.  Schematic  diagram  of  voltage  measurenent  circuit. 

FIGURE  2.5.  Measured  I(V)  curve  for  an  Al-PbBi  tunnel  junction  at 
1.815  K.  The  data  at  zero  voltage  is  supercurrent. 

FIGURE  2.6.  Electrical  model  of  the  tunnel  junction  including 
junction  capacitance  (C)  and  laser-induced  current  I^. 

FIGURE  2.7.  Observed  tangent  of  the  phase  shift  vs.  dV/dl  for  one 
tunnel  junction  at  three  different  temperatures.  The  straight-line 
fit  is  for  a  capacitance  of  2.4  X  10  F. 
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FIGURE  3.1.  Semiconductor-model  picture  of  disallowed  (3.1a)  and 
allowed  (3.1b)  tunneling  from  a  state  with  energy  for  two 
slightly  different  voltages  V. 

FIGURE  3.2.  G(E,V)  dE;  The  contribution  a  quasiparticle  in  the 
energy  range  E  to  E  +  d£  makes  to  the  tunneling  current  I(V). 
Parameters  A  and  A  2  taken  as  1.475  and  0.259  meV 

respectively. 

FIGURE  3.3.  Experimental  normalized  deconvolution  elements 
Dijlj/f^  vs.  V  for  one  particular  bin  "i".  The  temperature  is 
assumed  T  ■  1.586  K,  A  ”  1.475  meV,  and  A2  *  0.259  meV.  The 
energy  bin  is  [0.50,  0.55]  meV. 

FIGURE  3.4.  Unfolded  equilibrium  values  of  f(E)  vs.  E.  The  solid 
line  is  a  fit  to  the  Fermi  distribution  with  T  «  1.506  K.  Unfolding 
parameters  are  A  •  1.442  meV,  A2  *  0.2053  meV,  T  -  1.506  K, 
R  *  51.3JX  ,  and  g  ■  1/200  kA. 

FIGURE  3.5.  Unfolded  equilibrium  values  of  f(E)  vs.  E  normalized 
to  the  Fermi  distribution  with  temperature  1.506  K.  Parameters  are 
as  in  Fig.  3.4. 
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FIGURE  3.6.  Relative  shapes  of  I{V)/I(V)  predicted  for  changes  in 
Ai,  A2»  parameter  is  varied  per  curve.) 
Equilibrium  parameters  are  the  same  as  Fig.  3.4,  while 
I  5A]^  I  *  0.88  ^eV,  I  5A2  I  *0.88  |ieV,  |i  1.1  |ieV,  and 
St*  -  a  nK. 

FIGURE  3.7.  Unfolding  results  for  correct  values  of  parameters  A]^, 
A2»  (solid  dots),  and  with  errors  in  unfolding  parameters  (as 
marked) . 

FIGURE  3.8.  St*  model  unfolding  results  for  correct  parameters 
(dots)  and  with  errors  in  unfolding  parameters  (as  marked). 

FIGURE  3.9.  Semiconductor-model  of  tunneling  into  left-side  states 
with  energy  E^  for  two  different  voltages  V. 

FIGURE  3.10.  G(E,V);  The  contribution  quasiparticles  of  energy  E 
make  to  the  tunneling  current  I(V).  Parameters  are  A  and  E  are 
taken  as  0.55  meV  and  1.2  meV  respectively. 

FIGURE  3.11.  Normalized  deconvolution  elements  for  unfolding  I(V) 
for  a  symmetric  tunnel  junction.  The  temperature  is  assumed 
T  ■  2.17  K,  and  A  “0.55  meV.  The  energy  bln  is  {1.2,  1.25]  meV. 
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FIGURE  4.1.  Unfolded  nonequl librium  values  of  5f(H)  for 

aluminum-side  illumination.  Unfolding  parameters  are  meV, 

A 2*^ *2053  meV,  T»1.506  K,  R*51.3  ohms,  and  A 2“®*®®  |ieV. 

FIGURE  4.2.  Unfolded  nonequilibrium  values  of  ^f(E)  for 

aluminum-side  illumination  normalized  to  the  Fermi  distribution  with 

ic 

temperature  1.506  K.  The  solid  line  is  a  fit  to  a  T  model  with 
5X*“4nK.  The  dashed  line  is  a  |i*  fit  with  the  same  number  of 

•k 

excess  quasiparticles  as  the  T  fit. 

FIGURE  4.3.  Unfolded  nonequilibrium  value  of  aluminum  ff(E)  for 
lead-blsmuth-side  illumination  nonnalized  to  the  Fermi  distribution. 
The  equilibrium  parameters  are  essentially  the  same  as  Fig.  4.1, 
while  SA]^“  ^A2  *  -0.48  yeV.  Also  plotted  is  5f/f  expected  for 
a  T  distribution  with  oT  »  ZmK. 

FIGURE  5.1.  Semiconductor-model  picture  of  tunneling  between  two 
superconductors  at  zero  voltage.  Electron-like  (A)  and  hole-like 
(B)  tunneling  is  shown. 

FIGURE  6.1.  Detailed  electrical  schematic  of  the  SQUID  coupling 
circuit. 


FIGURE  6.2.  Tunnel  junction  layout  used  for  current  measurements. 
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FIGURE  6.3.  Schemaclc  dlagm  of  current  measurement  circuit. 

FIGURE  6.4.  Diagram  of  the  optical  path  from  the  laser  to  the 
optical  fiber. 

FIGURE  7.1.  Observed  laser-induced  voltage  shift  I  dV  I  vs.  V  for 
an  Al-PbBl  tunnel  junction  Illuminated  on  the  PbBl  side.  The 
temperature  was  1.588  K. 

FIGURE  7.2.  Measured  dV^  vs.  V  for  data  shown  in  Fig.  7.1. 

FIGURE  7.3.  Observed  dV*  vs.  dV/dl,  the  (voltage  dependent) 
dynamic  resistance  of  the  tunnel  junction.  The  straight  line  is  the 
best  fit  to  a  constant  thermoelectric  current  model. 

FIGURE  7.4.  Observed  I^  vs.  bath  temperature.  The  incident  laser 
power  was  held  fixed  at  IS  mW. 

FIGURE  7.5.  I^  vs.  incident  laser  power  for  4  fixed  temperatures, 
all  above  the  helium  lambda  point. 

FIGURE  7.6.  Fit  of  previous  data  to  theory.  Two  adjustable 

parameters  (c^^  and  71  )  were  chosen  to  fit  all  plotted  data. 
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FIGURE  7.7.  Plots  of  1  vs.  P  and  I  ^  vs.  P  for  an  Al~PbBi  txmnel 

o  o 

junction  at  2.29  K. 

LIST  OF  TABLES 

TABLE  A. I.  Description  of  terms  in  the  tunneling  hamlltonlan.  Both 


and  E^  are  defined  to  be  positive 


1 


I.  INTRODUCTION 

Considerable  experimental  and  theoretical  effort  has  been  spent 

studying  .the  nonequilibrium  quaslpartlcle  state  produced  by  an  external 

pair  breaking  mechanism  (Aronov  and  Splvak,  1978;  Langenberg,  1974; 

Langenberg,  1975).  Energy  externally  applied  to  a  superconducting  film, 

such  as  by  optical  Illumination,  perturbs  the  occupation  of  quaslpartlcle 

states  from  the  equlllbrltun  values  set  by  the  Fermi  function 

f(E)  ■  (1  +  Although  the  actual  distribution  may  be 

arbitrarily  complex.  It  Is  tempting  to  see  whether  It  can  be  described  by 

a  limited  number  of  easily  Interpreted  parameters.  In  an  early  attempt 

to  describe  nonequilibrium  states,  Owen  and  Scalaplno  (1972)  Introduced  a 

p*  model.  They  argued  that  emission  and  absorption  of  low  energy 

phonons  eventually  would  allow  the  quaslpartlcle  distribution  to  acquire 

the  temperature  of  the  phonon  bath,  while  a  slower  recombination  rate 

would  pose  a  bottleneck,  causing  a  buildup  In  number  of  quaslpartlcles 

which  could  be  taken  Into  account  by  a  shift  In  the  effective  chemical 

potential:  f(E)  ■  (1  +  e^^  I*  most  applicable 

for  weak  perturbation  at  low  temperatures’,  where  very  few  quaslpartlcles 

are  excited  and  hence  recombination  Is  slow.  Parker  (1975)  subsequently 

* 

showed  that  several  nonequilibrium  effects  could  be  explained  with  a  T 
model.  This  model  should  be  valid  for  those  cases  where  thermallzatlon 
by  the  low  energy  phonons  Is  the  slowest  step,  so  that  recombination  and 
pair  breaking  by  recombination  phonons  dominate  the  energy  relaxation. 
The  relatively  poor  phonon  coupling  to  the  low  energy  phonons  may  allow 
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Che  qua slpar tides  and  high  energy  phonons  to  reach  some  new  effective 
temperature  T*  significantly  different  from  the  bach  temperature.  In 
this  case  the  quasiparticle  occupations  could  be  approximated  by 

it 

E/ltT  —1 

f(E)  >(!■>■  e  )  .  Of  course,  one  can  combine  these  approaches, 

using  both  and  T*  to  give  a  2-parameCer  fit  to  Che  distribution 

function.  In  a  less  simplistic  approach  to  the  problem  Chang  and 
Scalaplno  (1978)  succeeded  in  solving  Che  linearized  kinetic  equations 
for  different  degrees  of  phonon-to-bath  coupling  and  several  special 
cases  of  excitation,  including  one  which  resembled  optical  illumination. 
For  poor  coupling  to  Che  bath,  they  found  their  results  closely  resembled 
T  "like  distributions. 

Quasiparticle  tunneling  has  been  used  in  much  of  Che  detailed 
examination  of  Che  validity  of  the  theoretical  models.  (See  Fig. 
l.I.)  Pioneer  work  by  Parker  and  Williams  (1972)  measured  the  gap 
reduction  under  optical  illumination,  which  they  Interpreted  in  terms  of 
p*.  Hu  et  al.  (1974)  were  able  to  actually  observe  the  dynamic 
relaxation  in  Che  number  of  quasiparticles  by  observing  Che  time 
dependence  of  the  quasiparticle  tunneling  current  after  pulsing  a  film 
with  a  laser.  They  succeeded  in  measuring  the  quasiparticle  relaxation 
time  for  tin  but  did  not  find  evidence  of  a  first-order  phase  transition 

it 

to  Che  normal  state  as  was  predicted  within  the  p  model. 

This  work  will  discuss  two  topics  related  to  illuminated 

superconducting  tunnel  Junctions.  A  new  method  for  using  asymmetric 


exponential  at  small  voltages. 
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tunnel  junctions  as  a  sensitive  and  accurate  probe  of  the  superconducting 
state  will  be  presented  and  will  be  applied  to  experifflental  results  for 
optically  illuminated  aluminum-aluminum  oxide-lead  bismuth  (Al-PbBl) 
tunnel  junctions.  Results  tend  to  confirm  the  validity  of  the  T  model 
to  describe  the  nonequilibrium  quasiparticle  occupation  function  In  the 
system. 

Also  presented  will  be  the  discovery  of  a  new  thermoelectric  effect 
In  tunnel  junctions.  Experimental  observations  confirm  the  existence  of 
a  thermoelectric  current  present  through  tunneling  barriers  between 
metals  at  different  temperatures.  The  existence  of  this  thermoelectric 
effect  can  explain  puzzling  discrepancies  between  experimental  results 
and  theoretical  predictions  In  several  earlier  superconducting 
thermopower  experiments. 


/  . 
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II.  EXPERIMENTAL  TECHNIQUES 


Accurate  unfolding  of  quaslpartlcle  I(V)  curves  to  obtain  occupation 
numbers  required  development  of  a  number  of  techniques  new  to  our 
laboratory.  First,  a  method  for  successfully  making  tunnel  junctions  had 
to  be  found  which  would  work  with  the  available  equipment.  The  tunnel 
junctions  had  to  display  very  sharp  gap  strtxctures  and  have  reasonable 
resistances.  Sections  A  and  B  describe  the  methods  I  used  to  fabricate 
high  quality  tunnel  junctions  reproducably .  Second,  the  I(V) 
characteristics  must  be  taken  with  sufficiently  high  accuracy  that  Input 
noise  would  not  make  the  output  occupation  numbers  meaningless.  It  was 
quickly  determined  that  taking  data  from  XY  plots  of  I(V)  was  both 
tedious  and  Insufficiently  accurate.  Sections  C  and  D  describe  a 
microcomputer-based  system  I  designed  and  built  to  control  and  monitor 
Che  laser  illumination  experiments.  Finally,  section  E  provides  details 
on  thermometry  techniques. 

A.  Substrate  Cleaning  Procedure 

For  the  results  reported  here,  polished  single  crystal  sapphire 
substrates  1  X  1/4  X  O.OIS”  In  size  were  obtained  from  Adolph  Heller  Co. 
The  substrates  were  each  recycled  many  times.  Careful  cleaning  of  the 
substrates  was  found  to  be  essential  for  making  uniform  metal  films.  The 
substrate  cleaning  procedure  given  below  evolved  from  a 
trial-and-many-errors  search  for  a  method  of  obtaining  consistent  film 
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quality. 


1.  Used  substrates  were  placed  In  a  beaker  of  distilled  water. 
Potassium  hydroxide  (KOH)  flakes  were  added  to  just  cover  the 
substrates.  The  KOH  dissolved  away  metal  films  and  loosened 
Indium  dots  from  earlier  sample  junctions.  Stirring  once  or 
twice  over  a  half  hour  period  usually  was  sufficient  to  clean 
off  the  substrates,  otherwise  this  step  was  repeated. 

2.  The  substrates  were  rinsed  thoroughly  with  distilled  water. 

3.  The  substrates  were  put  in  a  2  percent  solution  of  "Micro" 
cleaner  (manufactured  by  International  Products  Corp.)  which 
was  then  brought  to  a  boil  over  a  Bunsen  burner.  Heating  the 
solution  served  to  speed  up  the  degreasing  action.  The  hot 
mixture  was  then  ultrasonlcally  agitated  for  a  minute. 

4.  The  samples  were  again  thoroughly  rinsed  In  distilled 
water,  this  time  with  ultrasonic  agitation. 


5.  The  substrates  were  allowed  to  dry  In  air  under  a  Tensor 
lamp.  The  substrates  could  then  be  stored  until  needed. 
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B.  Tunnel  Junction  Fabrication 

Three  substrates  at  a  tloe  were  mounted  on  the  bottom  of  a  copper 
block  within  the  evaporator.  The  substrates  were  held  In  position  with 
small  dabs  of  Aplezon  M  high  vacuum  grease.  The  substrates  were  then 
covered  with  a  masking  plate  which  defined  an  exposure  pattern  for  the 
bottom  metal  layer.  The  evaporator  was  next  sealed  and  evacuated 
to  <  10  ^  Torr.  To  prevent  diffusion  pump  oil  backs treamlng  Into  the 
evaporator,  care  was  taken  to  be  sure  that  the  liquid  nitrogen  coldtrap 
was  always  filled. 

To  make  the  dirty  aliunlnum  layer,  It  was  necessary  to  Incorporate 
large  amounts  of  oxygen  Into  the  aluminum  during  evaporation.  Following 
a  suggestion  by  Kaplan,  evaporation  boats  of  Mo  overcoated  with  alumina 
(AI2O2)  were  used  to  hold  the  aluminum  pellets  to  be  evaporated.  As  the 
boat  was  heated.  It  outgassed  oxygen  In  addition  to  evaporating  aluminum 
from  the  pellets.  The  chamber  pressure,  as  measured  by  an  Ion  gauge 
above  the  coldtrap,  would  rise  to  3-9  X  10  ^  Torr  during  evaporation. 
T^'s  of  the  films  made  In  this  manner  ranged  from  1.6  to  2.3  K. 

The  next  step  In  making  the  tunnel  Junctions  was  to  grow  an  oxide 
barrier  layer  on  the  aluminum.  Barriers  formed  In  air  at  one  atmosphere 
pressure  formed  Junctions  with  resistance  too  high  to  be  useable 
(1  k/lmm  )•  To  obtain  lower  resistance  Junctions,  a  small  dead  space  In 
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the  air  release  line  was  backfilled  with  oxygen  gas  to  3  pslg  (980  Torr) . 
When  this  controlled  volume  of  gas  was  released  Into  the  chamber,  a 
pressure  of  0.1  Torr  resulted.  This  gas  was  subsequently  pumped  out 
after  a  period  of  15  seconds  to  2  minutes. 

The  masks  were  changed  without  breaking  vacuum  by  using  a  mechanical 
feedthrough  to  slide  the  masking  plate,  aligning  a  second  set  of  masks. 
After  evacuation  of  the  chamber  to  10  ^  Torr  again,  a  Pb^  95^^0  05 
was  evaporated  to  complete  the  tunnel  junctions.  Figure  2.1  shows  the 
metal  pattern  of  a  completed  junction.  Each  substrate  contained  three 
cross-strip  tunnel  junctions.  The  bottom  aluminum  strip  was  common  to 
all  three  junctions.  A  total  of  eight  contacts  were  made  at  the  edge  of 
the  substrate  to  be  able  to  make  four-terminal  I(V)  measurements  on  the 
tunnel  junctions. 

For  optical  Illumination  studies.  It  was  desirable  to  be  able  to 
Illuminate  the  tunnel  junction  from  both  sides  of  the  substrate.  The 
design  of  the  subtrate  holder  used  to  accomplish  this  Is  shown  In  Fig. 
2.2.  The  substrate  was  held  by  a  tight  fitting  recessed  area  in  the 
bottom  block  of  the  substrate  holder.  A  top  block  was  then  fastened  to 
the  bottom  block,  completing  the  holder  assembly.  Both  pieces  were  made 
of  phenolic.  Large  central  cutout  slots  permitted  optical  access. 
Electrical  contacts  were  made  by  compressing  Indium  dots  (actually 
cylinders  of  Indium  0.25  mm  X  1.  mm  dla.  made  by  cutting  Indium  wire) 
between  the  substrate  lead  pads  and  tinned  flat  posts  epoxled  Into  the 
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figure  2.2.  Substrate 
block. 


top  block  of  the  substrate  holder.  The  indium  dot  method  of  making 
electrical  contacts  worked  without  a  single  failure  at  low  temperatures. 
An  obvious  additional  advantage  of  this  method  over  soldering  is  that 
heating  of  the  substrate,  which  can  easily  damage  the  junction,  is 
avoided  altogether. 

C.  The  Data  Acquisition  System 

After  mastering  the  black  arts  of  tunnel  junction  fabrication,  it 
quickly  became  apparent  that  taking  data  would  be  no  simple  matter.  In 
order  to  achieve  the  highest  possible  sensitivity,  a  phase  synchronous 
(lock-in)  measurement  scheme  was  employed.  The  laser  beam  was  chopped  at 
25  or  337  Hz,  and  the  synchronous  signal  was  detected  and  recorded.  For 
constant  current  biasing  the  observed  shift  in  voltage  would  typically 
amount  to  only  a  few  pV.  It  was  necessary  to  accurately  record  I(V)  and 
dV(V).  The  dV(V)  measurements  had  to  be  converted  back  to  dI(V)  by 
dividing  by  the  dynamic  impedance  of  the  tunnel  junction.  Carrying  out 
all  of  these  calculations  for  the  hundred  or  so  data  points  needed  to  do 
inversion  turned  out  to  be  very  difficult  to  do  by  hand. 

To  Improve  upon  our  conventional  XY  recorder,  pen-and-paper  data 
recording,  digital  technology  was  called  upon.  The  main  objective  was  to 
build  a  data  acquisition  machine  which  could  record  data  with  high 
accuracy,  good  resolution,  and  very  low  noise.  The  system  had  to 
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simultaneously  monitor  and  record  voltage  signals  from  voltage 
amplifiers,  lock-in  amplifiers,  ammeters  and  thermometers.  In  addition, 
it  was  desirable  that  the  system  adapt  readily  to  modifications  in  the 
design  of  Che  experiments. 

After  studying  the  available  prepackaged  data  acquisition  systems, 
it  was  decided  that  the  most  cost  effective  and  versatile  data 
acquisition  system  could  be  built  by  designing  around  a  microcomputer. 

The  microcomputer  system  built  centers  around  the  Digital  Equipment 
Corporation  (DEC)  LSI-11/2  family  of  components.  The  microcomputer 
consisted  of  a  backplane,  CPU,  memory,  disk  system,  and  serial 
input/output,  shown  schematically  in  Fig  2.3.  The  backplane  is  the 
chassis  which  electrically  connects  the  various  individual  printed 
circuit  boards.  In  order  to  allow  maximum  room  for  additions  and 
adaptations,  an  especially  large  backplane  (DEC  DDVll-B)  was  modified  to 
conform  to  the  LSI-11/2  wiring  scheme.  The  backplane  has  room  for  18 
cards  connected  onto  the  central  data  bus,  plus  additional  space  for  9 
more  cards  which  can  be  wired  individually. 

The  microcomputer  is  controlled  by  a  central  processing  unit  (CPU) 
capable  of  manipulating  16  bit  long  “words"  of  data  or  instructions-  The 
CPU  card  handles  memory  access  and  manipulates  the  available  data  coming 
from  teletype  units  and  other  input/output  devices.  The  main  memory 
storage  is  supplied  on  a  28K  (IK  *  1024)  word  semiconductor  memory  card. 
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More  permanent  storage  of  data  and  programs  is  done  on  double  density 
floppy  magnetic  disks,  each  capable  of  holding  1/4  million  words  of  data 
or  program  Information*  The  Charles  River  dual  drive  disk  unit  consists 
of  two  separate  read/wrlte  stations,  plus  Its  own  dedicated 
microprocessor  which  maintains  an  orderly  file  structure  and  checks  for 
read  errors.  The  computer  can  edit,  compile  and  run  standard  FORTRAN 
language  programs,  as  well  as  use  Its  own,  somewhat  more  efficient  MACRO 
assembly  language.  The  versatility  and  ease  of  use  of  the  computer  Is 
best  evidenced  by  the  number  of  group  members  who  compete  for  time  on  the 
machine  to  take  data,  run  time-consuming  theoretical  calculations,  or 
edit  manuscripts. 

The  ability  for  the  microprocessor  to  read  analog  signals  from 
amplifiers,  DVMs,  lock-ins,  etc.  Is  provided  by  Incorporation  of  an 
Analog  Devices  RTI-1250S  analog-to-digital  (A/D)  input  board.  The 
computer  accesses  the  board  as  If  it  were  a  memory  location.  Under 
program  control,  the  A/D  board  selects  any  of  eight  input  lines,  samples 
the  voltage  on  that  line,  converts  the  -10  V  to  10  V  input  to  a  12  bit 
Integer,  and  presents  that  digital  information  to  the  CPU.  For  added 
resolution  for  small  signals,  a  programmable  gain  amplifier  (PGA),  with 
computer  selectable  gains  of  1,2,4,  or  8,  Is  included  in  the  A/D  circuit. 
The  PGA  can  be  used  to  Increase  the  sensitivity  of  the  A/D,  at  the  cost 
of  decreasing  the  input  voltage  range.  Although  the  maximum  accuracy  of 
the  converter  remains  fixed  at  2  (0.05  percent)  of  the  full  scale 


reading,  the  resolution  can  be  set  as  low  as  0.5  mV  by  using  the  PGA.  In 
order  to  Increase  the  sensitivity  still  further,  two  of  the  analog  Input 
lines  are  connected  to  special  low  drift,  high  Impedance  amplifiers 
(Analog  Devices  S22B)  which  are  wired  to  provide  a  switch-selectable  gain 
of  1,  10,  100,  or  1000.  All  of  the  analog  Inputs  are  differential-input, 
so  that  "ground"  connections  for  analog  lines  coming  to  the  computer  are 
actually  measured  as  separate  voltages  and  are  automatically  subtracted 
from  the  signal.  In  this  manner  potential  ground  loops  through  the 
computer  are  avoided  and  line  pickup  along  the  cables  Is  eliminated  as  a 
problem. 

Complementing  the  A/D  board,  an  Analog  Devices  RTI-1252 
digital -to-analog  (D/A)  board  provides  analog  output  from  the  computer. 
Each  of  three  memory  locations  can  be  used  to  set  a  corresponding  D/A 
converter  to  any  voltage  between  -10  V  and  10  V  with  -  5  mV  resolution. 
The  D/A  features  are  useful  In  controlling  repetitive  tasks  during  the 
course  of  the  experiment,  such  as  adjusting  the  current  through  the 
tunnel  junction. 

To  reduce  60  Hz  noise,  the  dc  power  for  the  computer  and  support 
circuitry  is  regulated  to  have  very  low  line  ripple  (measured  as  0.001  V 
which  Is  quite  small  compared  to  the  0.5  V  ripple  typical  of  DEC  power 
supplies).  The  power  supply  was  shielded  and  located  away  from  the 
sensitive  analog  circuits. 
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D.  Experimental  Setup 

The  overall  system  schematic  is  shovm  In  Fig.  2.4.  The  current  was 
set  by  applying  a  known  voltage  through  a  large  resistance  (10  kA  to 
10  MA  depending  on  the  junction  resistance  and  the  desired  bias  point). 
The  current  was  measured  by  sensing  the  voltage  across  the  standard 
resistor  with  a  high  Impedance  isolation  amplifier,  which  then  fed  Into 
one  of  the  computer  A/D  Inputs.  The  voltage  across  the  Illuminated 
tunnel  junction  was  amplified  with  a  PAR  113  preamplifier  and  then  sent 
to  another  A/D  input.  The  recorded  current  voltage  characteristics  for 
an  Al-*PbBl  tunnel  junction  Is  shown  plotted  In  Fig.  2.5.  The  large 
range  In  measured  currents  (which  reflects  the  exponential  voltage 
dependent  conductance)  was  obtained  by  changing  the  bias  resistor  and 
taking  the  curve  piecewise.  Resolution  of  the  current  measurements  Is 
better  than  one  nanoampere. 

Two  lock-ln  amplifiers  were  used  to  record  both  phases  of  the 
laser-induced  voltage  shifts  across  the  current  biased  tunnel  junction. 
Junction  capacitance  Introduced  some  phase  shift  and  attenuation  Into 
measurements  of  the  laser-induced  dV .  Fig  2.6  shows  a  simple  electrical 
model  for  the  tunnel  junction  and  measuring  circuitry.  The  effect  of 
laser  Illumination  Is  similar  to  that  of  an  external  current  generator 
which  would  add  a  current,  dl^(V),  to  the  Junction.  For  dynamic  junction 
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at  zero  voltage 
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resistance  R  (voltage  dependent),  the  observed  output  voltage  dV(V)  Is 
given  by 

dV  -  dl^  R  /  (1  +  iwRC)  (2.1) 

where  U)  is  Che  2  ir  times  the  chop  frequency  and  C  is  Che  Junction 
capacitance.  The  junction  capacitance  decreases  the  observed  dV  signal, 
especially  for  large  dynamic  resistances. 

Corrections  can  be  made  for  capacitive  losses  if  Che  capacitance  is 
known.  Inspecting  Eq.  (2.1)  it  can  be  seen  that  the  observed  phase 
shift  angle  is  given  by 

tan(  e  )  -  -uiRC  (2.2) 

Fig.  2.7  shows  Che  tangent  of  the  measured  phase  shift  as  a  function  of 

Che  dynamic  resistance  of  Che  tunnel  Junction  for  a  single  tunnel 

Junction  at  three  different  temperatures.  The  straight  line  fit  implies 
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a  Junction  capacitance  of  2. A  X  10  F.  The  fit  value  of  the  capacitance 
could  then  be  used  with  Che  dynamic  resistance  computed  from  the  I(V) 
curve  Co  give  a  complex  Impedance,  which  would  finally  give  Che  correct 
dI(V)  for  a  measured  dV(V)  using  Eq.  (2.1). 

During  my  sample  runs,  a  FORTRAN  program  handled  signal  sampling  and 
setting  of  parameters.  One  of  Che  computer's  D/A  outputs  was  used  to 
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fit  Is 


22 


sweep  out  the  I(V)  tunneling  curves  in  steps.  For  each  data  point,  the 
computer  set  the  current,  paused  1/2  second  to  allow  for  relaxation  of 
transients  In  the  lock-ln  amplifiers,  and  then  accumulated  and  averaged 
separately  measurements  of  the  current,  voltage,  and  dV  for  the  tunnel 
Junction  over  an  averaging  period  of  a  second.  For  each  current  value, 
two  data  points  were  taken  on  the  upsweep,  and  another  two  were  taken  on 
the  downsweep  In  order  to  check  for  reproducibility  (no  problem)  and 
hysteresis. 


E.  Thermometry 


To  monitor  temperatures  during  precooling  of  the  rig,  simple  silicon 
diodes  (type  1N4004)  proved  to  be  convenient  and  Inexpensive  (8  cents 
apiece)  thermometers.  Biased  In  the  forward  direction  at  10  |iA,  the 
voltage  drop  across  the  diodes  were  measured  with  a  digital  voltmeter 
(DVM) .  At  room  temperature  the  voltage  drop  was  typically  0.40  V.  The 
voltage  Increased  quite  linearly  with  decreasing  temperature,  at  least 
down  to  liquid  nitrogen  temperature.  At  liquid  nitrogen  temperature  the 
voltage  drop  was  typically  1.00  V.  The  diodes  cycled  to  liquid  helium 
temperatures  without  failures,  and  had  temperature-voltage 
characteristics  which  were  quite  reproducible. 

At  liquid  helium  temperatures  the  voltage  across  the  diodes  would 
rise  to  roughly  2.6V.  The  20  dissipated  by  the  diode  was  sufficient 
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to  prevent  effective  vapor  cooling  bjr  helium  vapors  alone*  During 
transfer,  only  when  the  diode  was  covered  by  liquid  helium  would  the 
voltage  on  the  diode  remain  steady.  Thus  it  was  a  simple  matter  to  use 
the  diodes  as  level  detectors  as  well  as  precooling  monitors.  An  added 
advantage  of  use  of  the  silicon  diodes  was  that  the  voltage  drops  across 
the  diodes  were  much  larger  chan  characteristic  thermal  emf's  or  ohmic 
voltage  drops  along  the  leads,  so  null  measurements  and  f our-termlnal 
measurements  were  both  unnecessary. 

To  measure  the  exact  temperature  during  data-taklng  runs,  later 
experiments  used  a  commercial  Cryocal  CR-IOOO  calibrated  doped  germanium 
resistance  thermometer.  Care  was  taken  to  have  the  thermometer  at  the 
same  height  as  the  tunnel  junctions  and  as  close  as  possible  to  Che 
samples.  The  bath  temperature  could  be  monitored  to  within  a  few  mK 
using  a  Kelthley  160  DVM  to  measure  the  voltage  across  the  constant 
current  biased  germanium  element.  The  DVM  analog  output  was  then  sampled 
by  the  computer,  which  could  convert  the  measured  voltage  into  a 
temperature  using  a  lookup  table  stored  in  memory. 


III.  UNFOLDING  PROCEDURE 


The  potential  usefulness  of  tunneling  data  to  determine 
nonequilibrium  occupations  was  demonstrated  by  Chang  and  Scalaplno 
(1976),  who  calculated  tunneling  curves  for  occupation  functions  which 
might  result  from  several  forms  of  perturbation.  Kaplan  et  a^.  (1977), 
Chang  (1977),  and  Klrtley  et  al.  (1978)  subsequently  studied  hole-llke 
and  electron-llke  Injection  of  quaslpartlcles  In  tunnel  junctions  and 
compared  the  observed  current  to  models  of  branch  Imbalance  to  extract 
relaxation  times.  In  a  series  of  experiments  using  a  double  tunnel 
junction  structure,  Gray  and  Wlllansen  (1978,  1979a,  1979b)  unfolded  I(V) 
data  to  observe  significantly  nonthermal  distributions  Induced  by  Intense 
quaslpartlcle  Injection.  Studying  pair  breaking  tunneling,  Balkashln  et 
al.  (1979)  found  non-T*  behavior  for  low  Intensity  illumination  of  lead 
tunnel  junctions. 

Most  recently  Jaworskl  and  Parker  (1979)  reported  results  from 
optically  Illuminated  symmetric  tin  tunnel  junctions,  using  a  matrix 
unfolding  technique  to  extract  occupation  numbers.  They  were 
unsuccessful  in  directly  calculating  the  occupation  function,  but  did 
perform  a  multi-parameter  fit  In  powers  of  I/E. 

This  chapter  presents  an  improved  method  (Smith  £t  ^. ,  1980)  for 
using  asymmetric  tunnel  junctions  as  a  <  ^nsltlve  and  accurate  probe  of 


the  superconducting  state. 
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A.  Derivation  of  Unfolding  Technique 

Figure  3.1  shows  Che  semiconductor  model  picture  of  quasiparticle 
tunneling  in  an  asymmetric  tunnel  junction  biased  at  voltage  V.  For  a 
particular  aluminum  quasiparticle  state  with  energy  the  voltage  has 
been  chosen  just  less  than  No  tunneling  from  this  state  is 

allowed  because  Che  quasiparticle  would  fall  into  Che  gap  of  Che  other 
superconductor.  (The  same  is  true  for  the  complementary  hole-like 
quaslparclcle  state  on  the  lower  branch  with  Che  same  energy.)  If  the 
voltage  is  increased  just  slightly,  as  shown  in  figure  3.1b,  tunneling 
will  be  allowed  for  the  upper  branch  states. 

Algebraically  Che  semiconductor  model  can  be  written: 

I(V)  -  (1/eR)  f  N^(E)  N2(E+eV)  [f^(E)  -  fjCE+eV)]  dE  (3.1) 

-«» 

(£  script  1  referring  to  the  lead  bismuth;  subscript  2  referring  to  the 
aluminum).  The  effect  of  the  gap  in  the  density  of  states  N^(E)  is  such 
that  aluminum  quasiparticles  of  a  given  energy  E  will  not  contribute  to 
the  observed  current  if  eV  <  Aj  ~  E*  The  sharp  onset  of  tunneling  from 
these  states  at  the  voltage  (Aj  ~  E  )/e  provides  the  means  to  recognize 
Che  extent  to  which  quaslparticle  states  of  energy  E  are  occupied. 
Analysis  becomes  more  complex  for  high  energies  (E  >  A^^)  because 
tunneling  from  these  states  does  not  have  the  same  sort  of  sharp  onset  as 
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is  Che  case  with  the  lower  energy  states*  Fortunately  these  highest 
energy  states  are  sparsely  populated  at  the  temperatures  of  interest  and 
therefore  contribute  little  to  observed  currents.  It  Is  convenient, 
therefore,  to  restrict  unfolding  analysis  to  calculation  of  aluminum 
quaslparclcle  occupations  of  states  with  energy  less  than 

The  tunneling  equation  may  be  manipulated  Into  a  more  convenient 
form  for  Inversion.  Since  lead  bismuth  states  with  E  <  A 
nonexistent,  it  Is  possible  to  set  fj^(E)  ■■  f2(E)  for  low  energies  and  to 
define  their  common  value  as  f(E).  Because  the  contribution  of  high 
energy  quasiparticles  to  the  current  is  insignificant  within  the  stated 
constraints,  it  Is  possible  to  extend  the  occupation  equality  requirement 
CO  Include  higher  energies  also,  with  no  noticeable  error.  Tunneling 
effects  due  to  branch  imbalance  (Gray  and  Wlllemsen,  1979)  should  be 
negligible  because  Che  optical  pair  breaking  creates  essentially  the  same 
number  of  electron-like  and  hole-like  excitations.  This  allows  use  of 
the  identity 

f(E)  +  f(-E)  -  1  (3.2) 


to  eliminate  Che  dependence  on  negative  energies.  Writing  equation  (3.1) 
in  terms  of  positive  energies,  and  regrouping,  yields: 


I 


28 


I(V)  -  (1/eR)  /  [N^(E)N2(E-eV)  -  N^(E-eV)N2(E)  + 

o 

Nj^(E+eV)N2(E)  -  Nj(E)N2(E+eV)]  f(E)  dE  (3.3) 

where  Nj^(E)  Is  defined  to  be  zero  for  E  <  Note  that  for  E  <  . 

and  eV  <  “  A2  »  only  the  third  term  is  nonzero.  (In  order  to  avoid 

pair-breaking  tunneling,  data  should  be  taken  only  for  V<  (Aj^ 

Because  biasing  problems  made  measurements  for  voltages  in  the  negative 
resistance  region  (A  “A2)/®  <  V  <  (A  ^  +A2)/®  impractical,  the 
tighter  restriction  that  V  <  (Aj^  was  actually  imposed  for  the 

data  used  for  unfolding.)  Following  the  notation  of  Javorski  and  Parker 
(1979),  the  convolution  function  G(E,V)  may  be  defined  such  that 


I(V) 


I 

/ 


G(E,V)  f(E)  dE 


(3.4) 


In  theory  the  linear  expression  for  I(V)  in  terms  of  f(E)  may  be 
Inverted  to  give  f(E)  from  observed  I(V).  The  objective  is  to  accomplish 
the  unfolding  with  high  accuracy  in  f(E)  with  sufficient  energy 
resolution  to  study  departures  from  equilibrium.  For  energies 
E  <  A]^  '  A2>  ^^®  inversion  process  is  well  behaved,  l.e.  I(V)  is  very 
sensitive  to  fractional  changes  in  f(E)  and  vice  versa.  Figure  3.2  shows 
the  behavior  of  the  convolution  function  G(E,V)  for  a  set  of  different 
values  of  E.  The  sharp  cutoff  for  eV  <  A^  “  E  suggests  that  I(V) 


mi 


1.0  (mV) 


FIGURE  3.2.  G(E,V)  dE:  The  contribution  a  quasiparticle  in  the 
energ7  range  E  to  E  +  dE  oiakes  to  the  tunneling  current  I(V). 
Parameters  and  A  2  taken  as  1.475  and  0.259  meV 

respectively. 
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measurements  for  voltages  near  (A ^  ”  E)/e  will  be  most  Important  in 
determining  the  occupation  number  £(£)> 

In  practice  it  is  convenient  to  remove  the  strongest  energy 
dependence  of  f(£)  by  incorporating  the  assumed  Fermi  function  into 

a  new  normalized  occupation  function  f  defined  \sy. 

f(E)  -  f(E)/f°(E)  (3.5) 

A  corresponding  G(E,V)  is  defined: 

G(E)  -  G(E)  f°(E)  (3.6) 

which  preserves  the  form  of  the  convolution 

r(V)  -  /  G(E,V)  f(E)  dE  (3.7) 

o 

The  inversion  process  is  made  tractable  by  taking  measurements 
(j“l»»-«.n)  and  dividing  the  energy  spectrum  into  "bins" 

[Ei.Ei^l],  (i«l,...,m).  For  a  given  energy  bin  f(E)  may  be  replaced  by 
a  mean  value  and  G(E,V)  may  be  integrated  to  produce  a  form  convenient 
for  applying  discrete  matrix  methods. 
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fi  *  raged  over 


^ji 


r"'  - 

J  G(E,V  )  dE 


(3.8) 


(3.9) 


Incorporating  matrix  notation  gives: 


■  I'ji  ■  (  i  i  >j 


(3.10) 


or  more  compactly: 


I  -  G  f 


(3.11) 


If  the  number  of  data  observations  n  exactly  equals  the  number  of 
unfolding  bins  m,  then  equation  (10)  represents  the  n  Independent  linear 
equations  which  can  be  solved  by  matrix  inversion  to  obtain  the  f 
values.  However  if  more  data  points  are  taken,  the  inversion  becomes 
overdetermined,  which  allows  the  advantageous  use  of  a  least  squares 
unfolding  procedure  with  resultant  noise  averaging.  For  the  unfoldings 
reported,  typically  45  data  points  were  used  to  calculate  20  f  values, 
so  the  unfolding  was  well  overdetermined. 


i. 
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To  solve  the  set  of  overdetermined  equations,  a  least  squares 
formula  was  developed.  The  mathematical  criterion  that  the  choice  of  the 
20  f  values  must  minimize 


-  X[(l(Vj)  -£G^^fp/  (3.12) 

J 

was  used  to  specif ly  the  unfolding.  The  meastirement  uncertainty  for  each 
data  point  a  has  been  Introduced  to  correctly  weight  the  data  points. 
For  given  Input  currents  and  uncertainties,  and  computed  C  values,  the 
function  describes  a  simple  paraboloid  In  20>space,  and  hence  has 
only  one  minimum.  Standard  procedures  (Bevlngton,  1969)  solve  for  the 
least  squares  fit  values: 

f  -  [(^  W  G)  ^  W]  I  -  D  I  (3.13) 

(1/  <Tj)^  (3.14) 

where  W  Is  a  diagonal  weighting  matrix  to  account  for  Instrument 
resolution  Oy  and  D  Is  the  resultant  deconvolution  matrix.  Of  course 
for  the  Inversion  not  to  be  undefined  the  number  of  data  points  n  must  at 
least  equal  the  number  of  bins  m. 
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A  coefficient  may  be  constructed  from  D  to  help  understand  which 
current  measurements  contribute  most  highly  In  determining  occupation  of 
a  given  energy  range.  The  normalized  elanent  Indicates  the 
fractional  contribution  of  the  current  measurement  to  the  calculation 
of  the  1^^  occupation  number.  Figure  3.3  shows  the  deconvolution 
fractions  /f^  for  an  experimental  equilibrium  set  of  data.  As 
expected,  the  current  measurements  which  are  most  significant  In 
determining  f(E)  are  those  taken  at  voltages  near  (A  ~  E)/e. 

The  accuracy  and  sensitivity  Inherent  In  the  unequal  gap  unfolding 
procedure  is  Illustrated  In  figures  3.4  and  3.5,  which  show  unfolding 
results  for  a  set  of  equilibrium  data.  The  parameters  Aj^,  A  2  T 
have  been  adjusted  to  obtain  the  best  fit  consistent  with  the  location  of 
the  A^  +A2  current  step;  it  was  not  found  to  be  necessary  to  include 
a  gap  smearing  parameter  to  describe  departures  from  the  BCS  density  of 
states.  However  It  did  slightly  improve  the  fit  to  correct  for  a  small 
parallel  constant  leakage  conductance  typically  4000  times  smaller  than 
that  of  the  normal  state.  Mean  deviation  relative  tc  the  Fermi 
distribution  Is  less  than  two  percent  over  three  orders  of  magnitude  In 
f(E),  while  the  absolute  resolution  for  equilibrium  quaslpartlcle 
occupation  Is  In  the  10  ^  range  for  the  higher  energy  states. 


/  , 


-0.5' 


1.0  (mV) 


FIGURE.  3.3.  Experimental  normalized  deconvolution  elements 
D. .I./f.  vs.  V  for  one  particular  bln  "1".  The  temperature  Is 
asiuied  T  -  1.586  K,  Ai  -  1-475  meV,  and  A2  “  0.259  meV.  The 
energy  bln  Is  [0.50,  0.55]  meV. 


paraneters  are  Ai  *  1<442  aeV,  A?  ~  0.2053  neV,  T  ■  1.506  K, 
R  -  51.3il  ,  and  g  -  1/200  kil. 
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FIGURE  3.5.  Unfolded  equilibrium  values  of  f(E)  vs.  E  normalized 
to  the  Fermi  distribution  with  temperature  1.506  K.  Parameters  are 
as  In  Fig.  3.4. 


B.  Unfolding  for  small  perturbations 


For  the  case  of  small  nonequilibrium  unfolding,  where  only 
departures  from  equilibrium  currents  are  measured,  a  small  modification 
of  technique  Is  necessary.  Irradiation  changes  f(E)  by  an  amount  f(E), 
which  Increases  the  observed  current.  An  Increase  In  f(E)  will  also 
have  the  effect  of  decreasing  the  gap,  which  will  subsequently  also  cause 
a  change  In  the  tunneling  current.  Differentiating  equation  (3.7)  yields 
the  linearized  approximation: 


SI(V) 


?0lE  (3.15) 


where  Is  the  shift  In  the  gap  parameter  for  superconductor  "1". 

Figure  3.6  shows  the  shape  of  the  incremental  change  In  current  resulting 
from  a  change  In  A  snd  a  change  in  A  2  compared  to  a  Sf(E)  of  the  |i* 
type  and  a  ^f(E)  of  the  T*  type.  These  shapes  are  best  understood  by 
reviewing  the  semiconductor  tunneling  picture  algebraically  described  by 
the  third  term  of  equation  (3.3).  Because  a  drop  In  A 2  decreases  the 
density  of  states  N2(E),  especially  near  the  gap,  the  effect  of  SA2 
$I(V)  Is  to  decrease  the  current,  partlculary  at  the  singularity.  A 
depression  in  A]^,  however,  at  a  given  voltage  V,  has  the  effect  of 
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FIGURE  3.6.  |elatlve  sjjapes  of  I(V)/I(V)  predicted  for  changes  in 
A  2^,  Aoi  |>  •  (Only  one  parameter  is  varied  per  curve.) 

Eq^llbrlua  parameters  are  the  same  as  Fig.  3.4,  while 
I  %Ail  -0.88  ^eV,  I  SAjI  -0.88  ^eV,  y  -1.1  yieV,  and 
St  *4  ^Ce 


r 


Initiating  tunneling  of  aluminum  quasiparticles  with  energies  In  the 
range  A  ^  ~  eV  -  Sa^  <  E  <A 1  ~eV.  Except  near  the  singularity,  a  drop 
in  A  contributes  fairly  uniformly  to  the  relative  Increase  in  current. 
Note  that  except  for  measurements  taken  at  voltages  near  A  ~^2 
very  difficult  to  differentiate  the  effects  of  a  drop  in  A  from  a 
change  in  the  effective  aluminum  chemical  potential  p*.  Consequently  it 
is  important  to  have  an  accurate  estimate  of  ^A^  in  order  to  make 
deductions  about  the  possible  p*  character  of  f(E). 

Assuming  the  equilibrium  state  f(E)  is  known,  it  is  possible  to 

correct  for  Sa  and  Sa  2  effects  by  subtracting  the  appropriate 

currents  contributed  by  the  respective  gap  drops.  The  resultant 
eff 

bI(V)  can  then  be  unfolded  in  a  similar  manner  as  for  I(V). 
Experimentally  the  quantity  SA^  ^A2  measured  by  observing 

the  current  jump  voltage.  It  is  necessary,  however,  to  make  some 
self-consistent  assumption  about  the  relative  sizes  of  S A  and  ^A2' 
As  a  check  on  the  Inferred  iA2»  resultant  Jf(E)  should  satisfy  the 
BCS  gap  equation  for  the  reduced  gap  valt;e.  Given  an  accurate 
parameterization  of  the  equilibrium  state,  current  represented  by  the 
iA^^  and  jA2  terms  in  equation  (3.15)  may  be  readily  subtracted  out  as 
a  background  current,  and  the  unfolding  may  proceed  as  before.  For 
comparison  to  theory,  it  is  often  then  useful  to  compute  Sf/f.  This  is 
readily  accompllahed  by  dividing  the  ac  Sf/f°  values  by  the  dc  f/f° 


results. 


AO 


The  advantages  of  using  the  ac  technique  are  that  the  high 
sensitivity  and  good  noise  rejection  of  phase  sensitive  detection  can  be 
used  to  measure  very  small  changes  in  occupation  numbers  (see  Chapter 
IV. )•  In  addition,  the  dc  (equilibrium)  data  can  be  used  to  precisely 
determine  the  Important  parameters  (Aj>A2>  needed  to  do  the 

unfolding. 


C.  Error  Analysis 

Before  applying  the  unfolding  technique  to  real  nonequilibrium  data, 
it  is  important  to  understand  what  types  of  errors  will  enter  into  the 
deconvolution,  and  what  the  sizes  of  these  errors  will  be. 

Systematic  unfolding  errors  can  result  from  incorrect  choice  of 
parameters  for  the  unfolding.  Although  location  of  the  gap  singularities 
provides  a  fairly  good  estimate  of  A  i  And  A2>  ^  known  from  normal 
state  resistance,  and  T  can  be  measured  with  a  thermometer,  small  errors 
in  parameter  estimates  can  still  present  problems.  For  the  unfolding 
procedure  to  provide  useful  information  about  the  nonequilibrium  state, 
it  is  important  chat  small  errors  in  input  parameters  should  not 
drastically  change  results  of  the  unfolding. 

To  check  the  sensitivity  of  the  output  f  values  to  erroneous 
choices  of  material  parameters,  a  theoretical  I(V)  curve  was  computed  by 
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numerical  Integration.  Figure  3.7  shows  the  results  of  unfolding 
calculated  BCS  model  I(V)  data  points  with  various  parameter  errors.  The 
dc  unfolding  is  sensitive  to  even  small  changes  in  the  input  parameters. 
This  sensitivity  can  prove  useful  in  getting  better  estimates  of  the  key 
parameters. 

For  the  nonequilibrium  unfolding,  the  results  are  highly  insensitive 
to  mistakes  in  the  input  parameters.  To  check  this  assertion,  Sl(V) 
values  were  calculated  using  a  model  T*  distribution.  The  nonequilibrium 
currents  were  then  unfolded  using  erroneously  selected  input  parameters. 
Figure  3.8  shows  the  results  of  unfolding  with  parameter  errors  similar 
to  those  used  for  figure  3.7.  The  excellent  error  immunity  is  the  result 
of  the  use  of  equilibrium  f  values  to  normalize  out  most  errors. 

A  related  source  of  systematic  error  is  miscalculation  of  the 
unfolding  elements  because  the  densities  of  states  for  the  two  metals  are 
not  perfectly  described  by  the  BCS  density  of  states.  In  reality,  there 
is  always  some  smearing  of  the  gap.  The  smearing  may  be  due  to 
Inhomogenelty ,  anisotropy,  or  finite  lifetime  effects  (Dynes,  e£  al. , 
1978).  For  both  PbBi  and  dirty  aluminum,  the  gap  width  is  quite  sharp. 
The  20  yeV  width  of  the  A2  ^Ise  does  provide  a  limit  to  the 
energy  resolution  obtainable  without  using  a  more  realistic  density  of 
states. 


E(meV) 

FIGURE  3.7.  Unfolding  resales  for  correct  values  of  parameters  Ai. 
Ao*  ®od  T  (solid  dots),  and  with  errors  In  unfolding  parameters  (as 
marked) . 
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E(meV) 

FIGOEIE  3.6.  St*  model  uafoldlng  results  for  correct  parameters 
(dots)  and  with  errors  In  unfolding  parameters  (as  marked). 
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Random  errors  enter  the  unfolding  from  three  sources:  electronic 
noise  In  measuring  the  current  of  a  data  point,  noise  in  measuring  the 
voltage  of  a  data  point,  and  roundoff  error  during  the  unfolding 
calculations. 

The  quantitative  effect  of  mlsmeasurement  of  current  can  be  readily 
calculated  using  equation  (3.12).  Assuming  independent  errors  in 
meas\irlng  the  current  Ij,  the  variance  in  unfolding  f^  is  given  by: 

^  J 

Because  constant  current  biasing  is  used,  the  most  important  random  error 
in  measuring  current  is  caused  by  the  limited  resolution  of  the 
analog-'to-dlgltal  converters  used  in  the  Analog  Devices  input  board  of 
the  microcomputer.  These  errors  have  been  measured  using  a  dummy  sample 
and  are  in  agreement  with  the  manufacturer's  specifications,  namely,  0.1 
per  cent  of  the  measured  current.  For  the  unfolding  Illustrated  in 
figure  3.3,  application  of  equation  (3.15)  suggests  an  uncertainty  in  the 
dc  f(E)  values  of  0.2  per  cent. 

Noise  currents  are  not  a  significant  problem  except  at  the  highest 
sensitivities.  Johnson  noise,  at  <  3  pA,  is  several  orders  of  magnitude 
smaller  than  the  externally  applied  current.  Ri  pickup  is  largely 
shunted  by  the  junction  capacitance,  while  line  pickup  is  minimized  by 
careful  shielding,  using  twisted  pairs,  and  avoidance  of  ground  loops. 
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Mlsmeasurement  of  the  voltage  of  a  given  data  point  can  be 
translated  Into  an  effective  error  in  reading  the  current  as 

-  (dl/dV)  a  (3.17) 

J 

The  voltage  error  is  primarily  due  to  low  frequency  ("1/f")  noise  in  the 
PAR  113  preamplifier  averaged  over  a  measurement  period,  and  amounts  to 
roughly  1  microvolt.  Because  the  current  is  roughly  exponential  in 
voltage, 

*  e  /kT  (3.18) 

For  T  ■  1.5  K,  this  error  is  0.75  per  cent*  Plugging  into  equation 
(3.16)  yields  an  error  in  f(E)  of  1.5  per  cent. 

The  final  source  of  random  noise,  computer  roundoff  error,  is  more 
difficult  to  calculate  directly.  The  individual  data  points  are  stored 
with  six  digit  accuracy  (far  more  accuracy  than  the  data  is  taken  with), 
but  the  iterative  arithmetic  operations  which  go  on  in  the  correlations 
and  matrix  inversion  each  introduce  roundoff  error.  To  test  the 
significance  of  these  errors,  computed  BCS  model  I(V)  data  points  were 
fed  in  as  simulated  data  and  unfolded.  The  results  were  accurate  to  the 
4  digit  print  out  resolution  used.  Roundoff  errors,  in  other  words, 
should  be  an  insignificant  component  of  the  total  noise. 


/ 
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D.  Unfolding  for  equal  gap  tunnel  junctions 

The  unfolding  procedure  described  In  sections  III.  A.  and  III.  B. 
Is  also  valid  for  symmetric  tunnel  junctions,  and  my  first  efforts  to  do 
unfolding  Involved  symmetric  tln-tln  oxlde-tln  tunnel  junctions. 
Unfortunately  there  are  significant  practical  limitations  to  useful 
application  of  these  techniques  (Jaworskl  and  Parker,  1979).  These 
Inherent  limitations  may  be  understood  within  the  framework  of  the 
notation  already  Introduced. 

For  purposes  of  discussion.  It  will  be  assumed  that  phonons  crossing 
the  oxide  maintain  f(E)  equal  on  both  sides  of  the  junction.  (Results 
vary  only  slightly  If  this  constraint  Is  eliminated.)  For  obvious  reasons 
the  restriction  to  study  only  those  states  with  E  <  ^  ^  will  be  dismissed 
for  this  symmetric  tunneling  case.  Furthermore,  the  necessity  of 
avoiding  pair-breaking  tunneling  Is  now  translated  to  the  requirement 
that  data  be  taken  at  voltages  eV  <  2A • 

Figure  3.9a  shows  the  semiconductor  model  picture  of  quaslpartlcle 
tunneling  In  a  s3nDmetrlcal  tunnel  junction  biased  at  a  voltage  V.  For  a 
particular  quaslpartlcle  state  with  energy  E^,  the  voltage  has  been 
chosen  just  less  than  (E^-^)/e.  The  electronrllke  quaslpartlcle  states 
noted  on  the  upper  branch  are,  by  the  zero  branch-imbalance  assumption, 
populated  the  same  as  the  hole-llke  quaslpartlcle  states  shown  on  the 
lower  branch.  As  shown,  both  electron-llke  and  hole-llke  tunneling  Is 
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allowed.  Because  the  lower  branch  Is  opposite  a  larger  density  of  states 
than  the  upper  branch,  the  net  current  flow  from  states  will  actually 
be  against  the  direction  of  the  applied  voltage!  The  situation  changes 
suddenly  as  the  voltage  Is  Increased  slightly,  as  shown  In  figure  3.9b. 
Tunneling  Into  the  lower  branch  states  Is  now  disallowed.  As  a  result, 
only  the  contribution  to  the  total  current  from  states  with  energy 
now  positive. 

Figure  3.10  shows  the  resultant  behavior  of  G(E,V)  for  this 
particular  value  of  the  energy  E^.  The  critical  voltage  Is  marked  as  V^, 

and  It  satisfies  the  condition  that  eV  -I*  E  *  A  •  There  Is  a  square  root 

singularity  as  the  gap  singularity  Is  approached  with  Increasing  voltage, 
followed  by  a  jump  to  a  fairly  steady  value  near  unity.  In  contrast  to 
the  asymmetric  tunneling  case  (see  figure  3.2)  where  tunneling  begins  at 
voltage  ( A 1  ~  E)/e,  symmetric  tunneling  changes  at  a  voltage  (E  -A)/e* 

The  major  difficulty  with  symmetric  junction  vinfoldlng  stems  from 
the  fact  that  the  current  jump  at  Is  masked  by  the  larger 

current  contributions  of  the  states  with  lower  energies  (and  hence  higher 
occupations!).  In  contrast,  with  asymmetric  timnellng  the  current 

Increase  for  a  given  energy  occurs  with  contributions  from  the  less 
populated  states,  which  can  be  readily  subtracted  out  (Gray  and 
Willemsen,  1979b). 


FICDRE  3.10.  G(E,V);  The  contribution  quesiparticles  of  energy  E 

make  to  the  tunneling  current  I(V).  Parameters  are  A  and  E  are 


taken  as  0.55  oeV  and  1.2  aeV  respectively. 
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Of  course  Che  mathematical  Inversion  process  is  still  correct  in 
theory.  For  noise-free  data  and  perfectly  defined  densities  of  states, 
Che  procedure  will  still  give  correct  f(E)  values  for  output.  For  real 
data,  however,  the  signal-to-noise  problems  become  unmanageable  (Jaworski 
and  Parker,  1979).  Quantitatively  this  can  be  understood  from 
calculation  of  the  unfolding  elements  I^/f^,  plotted  in  figure  3.11. 
The  large  values  and  slow  convergence  indicate  chat  Che  output  unfolding 
values  will  be  quite  sensitive  to  noise  in  Che  current  measurements. 
(See  equation  3.16). 

A  second  difficulty  with  symmetric  Junction  unfolding  involves 
handling  "leakage”  currents,  multiparticle  tunneling,  and  other  anomalous 
additional  currents.  High  sensitivity  is  needed  to  measure  Che 
occupation  of  the  high  energy  states.  Unfortunately  for  symmetric  tunnel 
junctions  Che  data  most  important  in  calculating  these  states  is  taken  at 
voltages  near  2^/e  where  these  anomalous  currents  are  largest.  For 
asymmetric  junctions,  on  Che  ocher  hand,  one  calculates  high  energy 
occupations  using  current  measurements  taken  near  V>0,  where  at  least 
some  of  Che  anomalous  currents  are  small. 

As  we  have  seen,  asymmetric  tunnel  junctions  offer  some  important 
advantages  over  symmetric  tunnel  junctions  for  unfolding  to  obtain 
occupation  numbers.  The  unfolding  procedure  can  be  used  to  extract 
highly  accurate  information  about  the  occupation  of  Che  nonequilibrium 


states 


FIGURE  3.11.  Nomalized  deconvolution  elements  for  unfolding  I(V) 


for  a  symmetric  tunnel  junction.  The  temperature  is  assumed 
T  *  2.17  K,  and  A  *  0.55  meV.  The  energy  bin  is  [1.2,  1.25]  meV. 
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Chapter  IV.  RESULTS  OF  NONEQUILIBRIUM  STATE  UNFOLDING 

A.  Aluminum-Side  Illumination 

The  unfolding  program  was  used  to  analyze  data  taken  for  both 
aluminum-side  and  lead-bismuth-side  illumination.  For  illumination  of 
the  dirty  aluminum  the  assumption  >>  ^^PbBi  appropriate. 

(For  clarity  the  subscripts  1  and  2  used  in  the  previous  chapter  have 
been  identified  explicitly  with  the  lead-bismuth  and  aluminum  films 
respectively.)  Photon  absorption  occurs  only  in  the  aluminum;  the 

lead-bismuth  is  perturbed  only  by  phonons  escaping  from  the  aluminum 

film.  Most  of  the  phonons  from  the  aluminum  film  have  energies  of 
roughly  (Chang  and  Scalapino,  1978).  The  lead-bismuth  film  is 

largely  transparent  to  these  phonons,  since  they  lack  sufficient  energy 
to  break  lead-bismuth  pairs.  The  lead-bismuth  film  should  remain  quite 
close  to  the  bath  temperature.  Therefore  I  assume  ^Ap^gj^  *  0  and 

determine  from  the  experimentally  measured  value  of 

£(Ap^g^  +  ^Al^*  assumption  will  be  tested  later.  The  values  of 

ApbBi>  A^2>  ^  values  fixed  by  fit  to  the  equilibrium  (dc) 

data. 

Results  of  unfolding  of  data  with  Illumination  on  the  aluminum  side 
are  shown  in  figures  4.1  and  4.2.  Increase  in  the  occupation  function 

was  limited  to  a  few  percent.  The  data  are  in  excellent  agreement  with 
Che  theoretical  jf(E)  resulting  from  a  rise  in  the  effective 


FIGURE  4.1.  Unfolded  nonequlllbrlum  values  of  6f(E)  for 
alumlmiffl-slde  llluiolaatlon.  Unfolding  parameters  are  meV, 

A2"0-2053  meV,  T-1.506  K,  R-51.3  ohms,  and  SA2"0*8® 


FIGURE  4.2.  Unfolded  nonequiXibrium  values  of  f(E)  for 
aluminum-side  Illumination  normalized  to  the  Fermi  distribution  with 
temnerature  1.506  K.  The  solid  line  is^a  fit  to  a  T  model  with 
5t  -4i*.  The  dashed  line  ^Is  a  y  fit  with  the  same  number  of 
excess  quasiparticles  as  the  T  fit. 
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quaslpartlcle  tanperature  by  an  amount  St*  equal  to  4  mK.  As  a 
self-consistency  check  on  the  unfolding,  the  fit  ST  value  can  be 
plugged  into  the  BCS  gap  equation,  i^ich  can  be  approximated  as: 

-  3.016  (1-t)  -  2.45  (1-t)^  (4.1) 

where  d  «  A/A(0)  and  t  *  T/T^.  (The  expansion  was  performed  by 

numerically  solving  the  BCS  gap  equation  and  fitting  the  results  for 
t<0.7  to  within  0.3  percent.)  For  the  data  plotted  here,  a  gap  drop  of 
0.88  +  0.05  peV  was  measured  from  shifts  in  the  (Ap^g^  +  A^2^^®  jump. 
This  compares  quite  well  to  the  computed  T  model  drop  of  0.84  peV. 

Thus  not  only  is  the  shape  of  f(E)  correct,  but  also  the  observed 

decrease  in  the  aluminum  gap  agrees  with  the  gap  decrease  expected  for  a 
T*  fit.  Furthermore,  the  earlier  assumption  that  8a  PbBi  “  °  appears 
justified.  Also  plotted  in  figure  4.2  is  Sf  expected  for  a  p* 

distribution  corresponding  to  the  observed  There  appears  to  be 

little  evidence  for  a  shift  in  chemical  potential  for  the  quasiparticle 
system. 

These  quantitative  unfolding  results  are  in  good  agreement  with 
conclusions  drawn  by  Gerbsteln,  et.  al.  (1978)  who  studied  intense 
optical  illumination  of  tin  in  a  tin-tin  oxide-lead  tunneling  structure. 
The  nonequilibrium  state  they  observed  closely  resembled  a  thermal  state 
of  higher  temperature  than  the  bath.  Several  potential  problems  involved 
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with  Intense  optical  illumination  are  avoided  with  weak  illumination. 
Inhomogeneous  state  effects  due  to  nonunlfotm  film  reflectivity, 
microscopic  laser  speckle,  or  material  Inhomogenelty  should  simply 
average  since  the  perturbation  from  equilibrium  Is  small  enough  to  be 
considered  In  the  linear  regime.  Also  the  relaxation  mechanism  and  rates 
should,  In  large  part,  be  equal  to  those  of  the  equilibrium  state  because 
the  excess  number  of  quaslpartlcles  Is  kept  small  compared  to  the  thermal 
occupation  values. 

In  order  to  compare  unfolding  results  to  theory.  It  Is  useful  to 
compare  Che  relative  speeds  of  Che  important  relaxation  channels.  Taking 
the  relaxation  times  of  Kaplan,  (1976),  the  recombination  race 

Is  roughly  three  times  faster  than  Che  thermal  phonon  scattering 
rate  Tg  One  does  not,  therefore,  expect  a  simple  model  to  give  a 
good  description  of  the  nonequilibrium  state.  The  phonon  crapping 
factor,  y,  (defined  as  the  inverse  of  the  probability  that  a 
recombination  phonon  will  escape  the  film  before  breaking  a  pair)  gives  a 
sense  of  what  fraction  of  the  recombination  phonons  break  pairs  and  hence 
are  able  to  adjust  Che  energy  distribution  of  Che  quaslpartlcles. 
Studies  by  Elsenmenger  et.  a^.  (1977)  show  the  transmissivity  of 

phonons  from  aluminum  to  sapphire  Is  only  0.14.  For  the  film  geometry 
used  here,  this  results  in  a  ^  of  roughly  20.  The  high  energy  phonons 
are  well  coupled  to  the  nonequilibrium  quaslpartlcles  as  is  assumed  under 
Parker's  T*  model.  Thus  for  our  experimental  conditions  the  T*  model 
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should  be  more  applicable  as  a  first  approximation  to  the  excited 
occupation  distribution,  in  agreement  with  out  experimental  results. 


B.  Lead-Blsmuth-Side  Illumination 


The  experiment  was  also  performed  with  illumination  of  the 

lead-bismuth  side  of  the  tunnel  junction.  From  the  perspective  of  the 
aluminum  film,  the  result  of  the  laser  illumination  is  to  create  a  phonon 
flow  from  the  lead-bismuth  film.  Numerical  studies  by  Chang  and 
Scalapino  (1978)  have  shown  that  most  of  the  energy  leaving  the 
lead-bismuth  film  is  in  the  form  of  recombination  phonons 

(■ft  Wpjji  ^^PbBl  *  ^  “eV).  This  energy  corresponds  to  roughly  ^2 
Although  the  energy  per  excitation  is  a  factor  of  1000  smaller  than  was 
the  case  with  the  direct  optical  Illumination  of  altimlnum,  the  general 
character  of  the  resultant  nonequilibrium  state  is  expected  to  be 
similar,  since  in  both  cases  the  excitation  energies  are  much  larger  than 

^Al' 

Because  of  the  relatively  poor  phonon  coupling  between  the 
superconducting  films,  an  assumption  that  »  S  Ap^g^  is  no  longer 

valid.  The  observed  shift  in  the  current  step  was 

again  used  to  gauge  the  value  of  ^Apgg^  +  Although 

quantitative  measurement  of  the  voltage  of  the  shifted  (Apgg^  ~  A^2^)/e 
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peak  In  current  was  impractical  due  to  hysteresis,  the  sign  of  the 
current  Increment  near  V  ~  observed  to  be  in  the 

opposite  direction  as  with  aluminum~slde  Illumination  (l.e. 

i  ^^PbBl  lead-bismuth-slde  Illumination).  The  value  of 
was  chosen  to  give  the  best  fit  to  a  T*-llke  nonequilibrium 
distribution  In  the  aluminum.  The  unfolding  results  are  shown  in  figure 
4.3  along  with  the  predicted  disequilibrium  for  a  St*  ^  2  mK 
distribution.  In  a  consistency  check,  the  alumlntim  gap  drop  predicted 
for  the  fit  ^T*  was  found  to  agree  with  the  parameter  Implied  by 

the  choice  of  ^Apgg^.  This  would  not  be  the  case  if  ^Ap^jgj^  were 
chosen  as  zero  as  was  the  case  for  aluminum-side  illumination.  It  Is 
tempting  to  attribute  the  small  deviations  from  the  T*  model  to  the 
effects  of  ^Apgg^  phonons  from  the  lead-bismuth;  however,  the  weakness 
of  the  self  consistency  check  on  the  parameters  for  this  case  is  such 
that  these  deviations  may  not  be  significant. 

The  jf(E)  distribution  shape  for  lead-bismuth  Illumination  Is 
consistent  with  results  for  aluminum  Illumination.  Although  the  Initial 
pair  breaking  energy  of  the  optical  photon  Is  much  greater  than  the 
the  large  phonon  trapping  factor  should  dominate  in 
determination  of  the  basic  nonequilibrium  jf(E)  shape. 

In  summary,  unfolding  of  I(V)  measurements  of  dissimilar  metal 
tunnel  Junctions  provides  a  useful  tool  for  measuring  the  quaslpartlcle 
energy  distribution.  Absolute  sensitivity  of  10  ^  and  relative  accuracy 
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FIGURE  4.3.  Unfolded  nonequilibrium  value  of  aluminum  £f(E)  for 
lead~blanuch-8lde  Ulumlnatlon  normalized  Co  Che  Fermi  dlacrlbuclon. 
The  equilibrium  paramecers  are  esaenclally  Che  same  as  Fig.  4.1, 
whl^e  SA^*  ^^2  *  Also  plocced  Is  (f/f  expecced  for 

a  T  dlscrlbuclon  wlch  (l  >  2mR. 
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of  a  few  percent  are  readily  obtainable  for  small  perturbations  from 
equilibrium.  Results  of  unfolding  experimental  data  from  optically 
Illuminated  alumlnum''~alumlnum  oxide~~lead-bismuth  tunnel  junctions  show 
the  aluminum  nonequl librium  state  is  well  described  by  a  T*  model 
behavior  over  several  orders  of  magnitude  in  £f(E).  This  last 
conclusion  is  valid  for  both  optical  and  high  energy  phonon  Illumination 
of  the  dirty  aluminum  film. 
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V,  THERMOELECTRIC  EFFECTS  IN  TUNNEL  JUNCTIONS 


A.  Introduction 


The  study  of  unfolding  leads  quite  naturally  to  the  stixdy  of 
thermoelectric  effects  in  tunnel  Junctions.  For  the  thin  oxide  barrier 
of  the  tunnel  junction,  temperature  gradients  of  10^  K/cm  could  be 
obtained  readily  with  laser  illumination. 

Thermoelectricity  in  superconductors  has  been  studied  in  several 
types  of  experiments.  A  temperature  gradient  across  a  bimetallic 
superconducting  loop  is  predicted  to  Induce  a  current  flow  around  the 
loop,  if  the  quasiparticles  have  nonzero  thermopower  S  (Gal'perin,  e£ 
al.  ,  1974;  Garland  and  Van  Harlingen,  1974).  A  number  of  different 
experimental  designs  have  been  used  to  look  for  such  an  effect,  most 
recently  a  laser-heated  loop  experiment  by  Schuller  and  Falco  (1978),  and 
a  toroidal  geometry  experiment  by  Van  Harlingen  and  Garland  (1978). 

A  second  group  of  experiments  has  concerned  thermoelectric  effects 
caused  by  interaction  of  temperature  gradients  and  supercurrents.  An 
effect  proportional  to  *^1  was  predicted  (Pethick  and  Smith,  1979)  to 
give  rise  to  a  local  charge  imbalance  within  the  quasiparticles  which 
could  be  measured  with  a  tunnel  Junction  probe.  The  magnitude  of  the 
effect  depends  on  the  mean  free  path  of  the  quasiparticles,  as  well  as 
A  >  but  is  entirely  Independent  of  the  normal  state  thermopower. 
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Experiments  by  Clarke  e£  al .  (1979),  Heldel  and  Garland  (1980),  and 

Van  Harlingen  (1980)  confirm  the  existence  of  an  emf  which  Is 

proportional  to  "v  *VT,  although  the  theoretical  magnitude  and 

s 

temperature  dependence  have  been  open  to  debate  (Schmid  and  Schon,  1979; 
Clarke  and  Tlnkham,  1980). 

This  chapter  and  following  chapters  will  treat  the  discovery  of  an 
entirely  different  type  of  thermoelectric  effect,  applicable  to  tunneling 
through  oxide  barriers  between  both  superconductors  and  normal  metals. 
The  effective  temperature  difference  across  the  oxide  barrier  gives  rise 
to  a  thermoelectric  current.  Observations  of  this  thermoelectric  effect 
will  be  presented  for  both  open  circuit  voltage  measurements  and  short 
circuit  current  measurements.  The  existence  of  this  thermoelectric 
effect  can  explain  discrepancies  between  experimental  results  and 
theoretical  predictions  In  several  earlier  superconducting  thermopower 
experiments  with  point  contacts. 

B.  Derivation  of  a  Tunneling  Thermoelectric  Current 


The  tunneling  barriers  In  tunnel  junctions  are  typically  10  -  40  S 
thick  oxides  of  one  of  the  electrodes.  Moderately  thick  oxide  barriers 
reduce  the  transmission  probability  of  electrons  to  cross  the  barrier  to 
-  lO”^^.  The  oxide  barrier  also  serves  as  a  barrier  to  phonon 
propagation  (Kaplan,  1979).  As  a  result,  the  metal  films  comprising  the 
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tunnel  junction  can  have  significantly  different  temperatures  in 
nonequilibrium  situations  (Melton,  et  al. ,  1980).  These  effects  are 
especially  important  at  low  temperatures,  when  the  phonon  heat 

4 

transmission  (  ~  T  )  is  low. 

Qualitatively  the  physical  processes  explaining  a  thermoelectric 
quasiparticle  current  at  zero  voltage  can  be  seen  in  Fig.  5.1,  in  which 
the  semiconductor-model  density  of  states  is  plotted  for  an  Al-PbBl 
ttmnel  junction.  A  heat  source  creates  an  equal  number  of  electron-like 
(upper  branch)  and  hole- like  (lower  branch)  quasiparticles  in  the  PbBi, 
so  that  the  superconductor  on  the  left  (PbBl)  is  effectively  at  a  higher 
temperature  than  the  superconductor  on  the  right  (Al).  The  Fermi 
energies  are  held  eqiul  externally  by  a  superconducting  shorting  wire. 
Electrons  tunnel  across  the  barrier  from  the  PbBl  to  the  Al  in  the 
tunneling  channel  labeled  A.  There  is  a  back  current,  labeled  B,  for 
states  in  the  lower  branch.  For  an  oxide  transmission  probability  X(E) 
that  is  energy  independent,  the  charge  transfer  for  the  two  processes 
exactly  cancels.  Harrison  (1961)  has  noted  that  the  densities  of  states 
of  the  metals  are  proportional  to  (dE/dk)  ^  and  the  group  velocities  for 
electrons  to  approach  the  barrier  are  proportional  to  dE/dk,  so  that  the 
branch  cancellation  is  independent  of  density  of  states  changes.  If 
transmission  probabilities  for  A  and  B  processes  are  not  exactly  equal, 
however,  there  will  be  a  net  thermoelectric  current,  which  can  be 
experimentally  measured  by  using  a  galvanometer  in  place  of  the 
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N(E) 


RE  5.1.  Semiconductor-model  picture  of  tunneling  between  two 
rconductors  at  zero  voltage.  Electronrlike  (A)  and  hole-like 
tunneling  la  ahown. 
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superconducting  shorting  wire* 

An  estimate  of  the  energy  dependence  of  the  electron  transmission 
probability  can  be  obtained  using  the  oner* band  WKB  approximation*  For  an 
oxide  barrier  potential  electron  energy  E  (both  measured  from 
and  oxide  thickness  s,  X(E)  is  given  by  (Hartman,  1964) 

ln[X(E)]  »  -cfs(«-E)’-^^  +  constant  (5.1) 


where 

OC  -  4ir(2m)^^^  /h  -  1*025  eV^^^  X’’’  (5.2) 

The  dependence  of  X  on  the  angle  of  Incidence  upon  the  barrier  has  been 
Integrated  out. 

We  will  restrict  our  analysis  to  study  of  tunneling  at  low  voltages, 
l.e.  V  «  ^/e*  In  this  case,  only  states  close  to  Ep  contribute 
significantly  to  the  tunneling  current*  The  transmission  probability  for 
these  states  may  be  expressed  in  terms  of  a  Taylor's  expansion  about  E^.: 

ln(X(E))  -  -b  +  Cj^  E  +  *  **  (5.3) 

The  fractional  change  In  X  as  E  varies  is  given  by 
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(1/X)  (dX/dE)  Ip  -  c,  -  d/dE  ln(X)  |  p  (5.4) 
F  F 


For  Che  model  described  by  Eq.  (5.1)  this  may  be  evaluated  as: 


-  0Cs/(2  ■  (23  +  ln(RA)]/(2^)  (5.5) 

2 

where  R  Is  In  ohms  and  A  is  In  cm  • 

Typical  values  of  Che  parameters  for  an  aluminum  oxide  tunneling 

barrier  would  be  s*10  &  and  0  *  3.7  eV.  (Konkin  and  Adler,  1974), 

Evaluating  equation  (5.5)  for  electron  states  near  Che  Fermi  surface  and 
2 

a  1  ohm  cm  junction  resistivity,  this  simple  model  gives: 


Cj^  •  3.1  eV 


(5.6) 


(It  should  be  noted  chat  Gundlach  (1973)  has  done  calculations  using  a 
slightly  more  sophisticated,  two  band  model.  This  model  Introduces  the 
effects  of  Che  valance  band  of  the  oxide  on  the  tunneling  probabilities. 
In  general  Che  results  of  these  calculations  are  to  lower  the  estimates 
of  C|^,  although  accurate  predictions  for  c^  are  rather  difficult.  The 
two  band  model  has  proven  superior  to  Che  one  band  model  when  applied  to 
explain  high  voltage  asymmetries  In  tunneling  I(V)  characteristics 
(Basavalah,  et  al.,  1974;  Ekrut  and  Hahn,  1980).) 
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If  the  relevent  energy  separation  of  the  two  branches  is 
(*  3  meV),  equation  (4)  suggests  that  the  A  currents  Illustrated  In 
figure  5.1  could  be  expected  to  be  on  the  order  of  one  percent  larger 
than  the  B  currents. 


The  exact  aagnltude  of  the  current  depends  on  the  excess  number 
density  of  excited  quaslpartlcles  that  are  able  to  tunnel  as  well  as 
details  of  the  junction  Itself.  The  semiconductor  model  predicts 


-  (l/eRX(0))  /X(E)  Nj(E)  N2(E)  [f^(E)-f2(E)  ]  dE  (5.7) 


where  is  the  predicted  zero  voltage  current,  Is  the  superconducting 
density  of  states  normalized  to  the  electronic  density  of  states  per 
spin,  N(0),  and  e  Is  negative.  The  normalized  transmission  probability 
can  be  expanded  as: 


X(E)/X(0)  -  1  +  Cj^  E  + 


(5.8) 


Using  this  expression  In  Eq.  (5.7),  the  first  term  Integrates  to  zero, 
leaving : 


Iq  -  (1/eR)  c^  /  E  N^(E)  N2(E)  [fj(E)  -  f2(E)] 


(5.9) 
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Although  the  semiconductor  model  used  to  describe  the  phenomena  of  a 
thermoelectric  current  correctly  predicts  the  rough  magnitude  and  correct 
sign  of  the  effect,  It  does  not  carefully  cake  Into  account  the  combined 
hole-electron  nature  of  quaslpartlcles.  (For  a  more  detailed  and 
rigorous  calculation  of  the  thermoelectric  current,  see  Appendix  II>)  The 
major  conclusion  should  be  valid,  however,  namely  that  a  thermoelectric 
current  proportional  to  Che  excess  number  of  quaslpartlcles  able  to 
tunnel  should  be  excited  across  the  oxide  barrier  of  the  tunnel  junction, 
i.e. 


I^cCn-n^(T)  (5.10) 

where  n(T)  Is  the  number  density  of  quaslpartlcles  (conventionally 
normalized  to  4N(0)A(0)I1»  XI  is  the  electrode  effective  volume,  and  n^ 
Is  the  thermal  equilibrium  value  defined  by  the  counter  electrode  (Al) 
effective  temperature.  (See  appendix  equations  (A. 7)  and  (A.8)). 

It  should  be  noted  that  the  thermoelectric  current  derived  here  is 
entirely  a  tunneling  barrier  effect.  The  current  is  independent  of  those 
electrode  material  parameters  normally  associated  with  thermoelectricity: 
the  thermopower  S,  the  quaslpartlcle  diffusion  length  .Aq,  the  mean  free 
path  etc.  The  only  Important  material  parameters  are  chose 
determining  Che  number  of  quaslpartlcles  which  are  able  to  tunnel. 
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C.  Dependence;  on  Bath  Temperature  and  Laser  Power 


Measurements  which  will  be  described  more  fully  In  Chapter  VII. 
were  taken  to  study  the  power  and  temperature  dependence  of  the 
thermoelectric  effect.  As  before,  a  laser  was  used  to  supply  heat  to  the 
lead-bismuth  film  of  the  tunnel  junction.  An  advantage  of  this  technique 
over  joule  heating  techniques  Is  that  It  avoids  Inductive  or  capacitive 
coupling  between  external  heaters  and  the  tunnel  junction.  A  difficulty 
of  the  laser  heating  technique  Is  that  the  effective  temperature  of  the 

•k 

Illuminated  film,  T  ,  must  be  estimated  rather  than  measured  directly. 

The  value  of  (n  -  n^)  Is  related  to  the  Incident  laser  power  P,  as 
well  as  the  gap  parameter.  The  number  of  quasiparticles  excited  may  be 
estimated  using  the  coupled  Rothwarf-Taylor  (1967)  rate  equations.  The 
steady  state  number  of  excited  quasiparticles  (In  normalized  units)  Is: 

n  -  -  (P  <E>  A  (0)  N(O)il)  (5.13) 

where  “^eff  Is  the  effective  lifetime  of  the  excitations  before 
recombination  occurs  and  <£>  Is  the  average  energy  per  quaslpartlcle 
••rker,  1975).  The  lifetime  is  itself  dependent  on  the  number  of 

i««. particles-  For  a  quaslpartlcle  to  recombine,  it  must  pair  with 
;uailparticle.  The  recombination  rate  for  quasiparticles  Is 


simply  proportional  to  the  number  of  combinations  of  pairs  of 
2 

quasiparticles  (n  /2),  so  the  effective  recombination  time  is  given  by: 

Tgff  “  n  (dn/dt)  ^  •  b  TT/n  (5.14) 

where  is  a  characteristic  time  for  the  superconductor  (Kaplan,  et 
al.  ,  1976),  b  is  a  dimensionless  constant  of  proportionality  (equal  to 

0.048),  and  "T  is  the  phonon  trapping  factor  (Kaplan,  1979),  which  may  be 
on  the  order  of  50  depending  on  Che  materials  at  the  interface. 

If  n  <<  1,  A  Is  Independent  of  n,  and  equations  (5.13)  and  (5.14) 
may  be  solved  self-consisCently ,  yielding: 

n  -  n^  -  {[  n^  +  P^b  r /«£>  A  N(0)yi  )]^^^  -  n^}/2  (5.15) 

Several  limiting  cases  of  Eq .  (5.15)  are  worthy  of  note.  Combining 

Eqs.  (5.15)  and  (A. 12)  gives  for  low  incident  power: 

«  (2/eR)  c^  P  T^'yb/(4n^N(0)il)  (5.16) 

The  current  is  simply  proportional  to  the  incident  laser  power,  with  a 
coefficient  chat  depends  on  temperature.  For  temperature-independent  c^ 
and  'd' ,  this  temperature  dependence  goes  as  n^(T) 
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At  low  temperatures,  where  n^(T)  becomes  small  exponentially,  the 
number  of  laser-excited  quaslpartlcles  may  exceed  the  number  of  thermally 
excited  quaslpartlcles.  For  sufficiently  large  laser  power,  Eqs.  (5.13) 
and  (A. 12)  lead  to 

a  (1/eR)  Cj^A  (P7T;/N(0)n)^''^  (5-17) 

The  linear  thermoelectric  current  saturates,  producing  a  square  root 
dependence  on  applied  laser  power. 

In  summary,  the  expected  power  dependence  predicted  by  this  simple 

Rothwarf-Taylor  equation  model  should  be  linear  for  small  laser  power  and 

moderately  large  temperatures.  For  sufficiently  low  temperature  (and 

therefore  small  n^) ,  the  predicted  thermoelectric  current  should  vary  as 
1  /2 

F  .  For  fixed  (small)  laser  power,  the  temperature  dependence  of  the 
current  should  roughly  follow  n^(T)  which  Increases  with  decreasing 
temperature.  All  of  these  predictions  are  based  on  temperature 
Independent  values  for  K  and  c^^. 
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Chapter  VI.  THERMOELECTRIC  EFFECT  MEASUREMENT  TECHNIQUES 

A.  SQUID  Circuit 

Measurement  of  the  thermoelectric  current  presented  challenging 
problems  in  low  noise  metrology.  To  measure  signals  down  to  1  pA,  at 
essentially  zero  voltage,  a  circuit  was  designed  around  a  S.H.E.  SQUID 
galvanometer.  The  quoted  current  noise  for  the  SQUID  system  was  20 
pA/(Hz)  .  In  order  to  improve  upon  this  noise  figure,  a 
superconducting  transformer  was  wound  using  superconducting  NbTi  wire. 
The  transformer  was  surrounded  by  a  lead  shield  to  reduce  magnetic  pickup 
from  external  fields.  The  in-circuit  current  gain  of  the  transformer  was 
measured  as  lA.  Use  of  the  computer  to  average  data  over  11  second 
Intervals  reduced  the  bandwidth  to  roughly  0.1  Hz.  The  total  system 
noise,  measured  by  taking  data  on  a  superconducting  tunnel  junction,  was 
0.6  pA,  which  is  very  close  to  the  theoretical  limit  of  O.A  pA  set  by  the 
SQUID-trans former  circuit. 

The  electronic  schematic  for  the  SQUID  circuit  is  shown  in  Fig.  6.1. 
The  SQUID  noise  bandwidth  was  limited  to  about  10^  Hz  by  including  a 
3  ohm  parallel  resistor,  Rp.  This  acted  as  a  low-pass  filter, 
attenuating  all  higher  frequency  signals.  The  actual  laser  signal  was 
chopped  at  25  or  337  Hz,  so  the  high  frequency  limit  presented  no 
difficulty.  Also  included  in  the  circuit  was  a  6  x  10  ^  ohm  series 
resistor,  Rg,  which  provided  a  low  frequency  rolloff  at  3  Hz.  This  was 
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done  CO  eliminate  crapped  flux  from  Che  SQUID  loop.  Note  Chat  in 
contrast  Co  Che  capacitive  corrections  required  for  Che  high  impedance 
experiments  described  earlier,  tunnel  junction  capacitance  has  no  effect 
on  these  low  impedance  measurements.  The  SQUID  circuit  effectively 
shorts  out  Che  capacitance. 

The  sample  geometry  was  modified  slightly  from  the  earlier 

experiments,  and  is  shown  in  Fig.  6.2.  The  tunnel  junction  (located  on 

2 

Che  upper  end)  was  somewhat  larger,  measuring  about  5  mm  .  Electrical 
contacts  to  Che  sample  were  made  by  Che  same  type  of  pressure  contact 
mount  as  was  described  in  chapter  2.  At  Che  ocher  end  of  the  sapphire 
substrate,  electrically  isolated  from  the  sample  junction,  a  simple 
aluminum  strip  was  evaporated  at  the  same  time  that  the  aluminum  was 
evaporated  for  the  tunnel  junction.  Four-terminal  measuranents  of  the 
isolated  strip  were  used  to  determine  T^  of  the  aluminum  electrode  of  the 
tunnel  junction  without  introducing  noise  into  Che  junction  circuit. 

The  overall  current  measurement  experiment  is  diagramed 
schematically  in  Fig.  6.3.  A  chopped  laser  beam  illuminated  the  tunnel 
junction,  generally  from  Che  lead-bismuth  side.  The  tunnel  junction  was 
coupled  through  the  superconducting  transformer  Co  Che  SQUID  circuit. 
The  SQUID  output  signal  was  filtered  with  an  Ithaco  electronic  filter  to 
remove  occasional  flux-jump  spikes.  The  output  from  the  filter  was  then 
fed  into  Che  two  lock-in  setup  described  earlier.  The  computer  was 
programed  to  sample,  average,  and  record  the  lock-in  output  signal  in 
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FIGURE  6.3.  Schematic  diagram  of  current  measurement  circuit. 
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addition  to  the  bath  temperature,  laser  po/er,  and  resistance  of  the 
aluminum  strip.  Runs  were  taken  varying  either  the  temperature  or  laser 
power. 


B.  Fiber  Optics 


The  extremely  high  sensitivity  of  the  SQUID  system  made  careful 
electronic  shielding  of  the  tunnel  junction  from  outside  rf  noise 
imperative.  Unfortunately,  a  major  source  of  the  rf  noise  in  the 
laboratory  was  the  electric  discharge  in  the  laser  tube.  Successful 
simultaneous  operation  of  the  SQUID  and  the  laser  was  found  to  require 
keeping  the  tunnel  junction  and  SQUID  in  an  rf  shielded  room,  while 
somehow  getting  the  laser  beam  into  the  shielded  room  and  down  onto  the 
tunnel  junction  sample. 

These  difficulties  were  overcome  through  the  use  of  an  optical 
fiber.  Recent  advances  in  optical  fiber  technology  have  produced  low 
loss  (  <10  dB/km) ,  low  cost  (^$I/meter),  commercially  available  fiber 
capable  of  transporting  powers  of  up  to  kilowatts  in  a  flexible  light 
pipe.  For  my  requirements,  I  wanted  to  be  able  to  transport  the  laser 
beam  over  a  distance  of  roughly  30  meters,  without  losing  more  than  a 
factor  of  ten  in  optical  power.  The  fiber  also  had  to  be  able  work  at 
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temperatures  as  low  as  1.5  K,  without  loss  of  mechanical  strength  or 
flexibility. 

Several  types  of  optical  fibers  were  considered  for  use  In  the 
experiment.  Plastic  fibers  were  ruled  out  because  of  their  high  loss 
(  >  1  dB/meter).  Fibers  of  silica  with  silicone  cladding  were  unsuitable 
because  the  silicone  changes  its  Index  of  refraction  at  low  temperatures, 
causing  the  fiber  to  become  quite  lossy.  A  difficulty  presented  by  fiber 
bundles  and  composite  fibers  was  differential  contraction  of  the  cladding 
material.  At  low  temperatures  the  differential  contraction  causes 
microbending  which  leads  to  high  optical  losses  Into  the  cladding 
material. 

All  the  above  difficulties  were  avoided  by  using  an  all  glass  fiber 
made  by  Coming  Glassworks.  Coming  was  generous  enough  to  respond  to  my 
request  for  a  small  sample  by  sending  5  km  of  their  Corguide  optical 
waveguide.  A  length  of  the  fiber  was  tested  In  liquid  nitrogen,  and 
found  to  be  strong,  flexible  and  almost  loss-free. 

There  were  two  major  constraints  to  a  practical  coupling  scheme  from 
the  laser  to  the  fiber.  The  first  involved  the  size  of  the  fiber.  The 
optically  conducting  core  of  the  fiber  was  only  63  ym  In  diameter.  The 
nearly  parallel  laser  light  could  be  focussed  down  to  a  spot  size  of 
f  X  ,  where  f  Is  the  f  number  of  the  lens  system,  and  X  Is  the  optical 
wavelength  (0.5145  ym) .  This  Imposed  a  constraint  that 
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f  <  100  (6.1) 

It  should  be  understood  that  this  limit  is  for  a  very  high  quality  lens. 
For  the  more  economical  lenses  actually  used,  this  limit  should  be  scaled 
down  somewhat.  In  general,  the  larger  f  is,  the  larger  the  spot  size  on 
the  fiber. 

The  second  constraint  on  the  laser-fiber  coupling  scheme  was  Imposed 
by  the  limited  acceptance  angle  of  the  fiber.  The  index  of  refraction  of 
the  fiber  material  falls  off  parabollcally  away  from  the  center  of  the 
channel.  This  serves  to  refocus  light  which  is  not  too  greatly  off  axis. 
For  light  inside  the  fiber  core  to  be  totally  Internally  reflected  by 
this  refocusing  mechanism,  the  departure  from  alignment  should  not  exceed 
the  critical  angle  set  by  the  index  of  refraction  profile  of  the  film. 
For  the  Coming  Corgulde  fiber  used  in  this  experiment,  it  was  necessary 
that  light  hitting  the  center  of  the  fiber  enter  from  an  angle  of  no  more 
than  0.2  radians  (  12°).  (For  light  near  the  edge  of  the  conducting  core 
the  critical  angle  is  even  smaller.)  This  imposes  a  constraint  that 

f  I  4.7  (6.2) 

A  tradeoff  was  made  between  small  spot  size  and  small  entrance  angle  by 
using  a  9  mm  focal  length  lens.  For  the  1.5  mm  diameter  laser  beam,  this 


80 


resulted  in  £/6.  Experimentally  this  was  found  to  couple  the  power  into 
the  fiber  most  efficiently,  for  the  argon-ion- laser  beam.  A  spot  size  of 
10  |jm  could  be  readily  obtained  and  scanned  across  the  end  of  the  fiber 
using  an  NRG  FB-1  ]CYZ  positioning  mount  for  the  fiber. 

A  practical  difficulty  in  the  use  of  the  fiber  was  that  its 
transparency  and  small  size  made  it  almost  invisible.  Although  the 
opClc2d  fiber  was  strong  enough  to  support  its  own  weight,  as  well  as 
take  some  fairly  rough  handling,  it  consistently  failed  when  it  was 
stepped  upon.  This  necessitated  running  the  fiber  as  close  to  Che 
ceiling  as  possible,  and  hence  out  of  harm's  way. 

The  fiber  was  run  from  the  XYZ  mount  up  to  the  ceiling,  through  a 
hole  in  the  wall,  across  another  room,  through  the  copper  screening  of 
the  rf  shielded  room,  and  down  to  Che  cryostat.  At  the  cryostat  it  was 
threaded  down  a  chin  teflon  Cube  which  was  then  sealed  at  Che  cop  with 
putty-like  Aplezon  Q  leak  sealant.  The  teflon  tube  served  to  protect  the 
fiber  from  abrasion.  The  fiber  was  positioned  above  the  tunnel  junction 
by  teflon-taping  it  onto  a  small  support  rod  located  on  the  substrate 
holder. 

The  power  transmission  for  Che  entire  system,  including  30  meters  of 
fiber,  was  measured  using  a  thermopile  detector  at  the  output  end  of  the 
fiber.  For  Che  optimized  system,  the  total  coupling  efficiency  (defined 
as  the  ratio  of  output  from  the  fiber  divided  by  Che  laser  ouput  power) 


was  0.65.  This  compares  quite  favorably  with  the  theoretical  limit  of 
0.78  Imposed  by  partial  reflections  (4  percent  per  surface)  at  the  glass 
surfaces  in  the  optical  path  (the  beam  splitter,  the  focusing  lens  and 
the  optical  fiber  itself). 

An  added  advantage  in  use  of  the  glass  optical  fiber  over  more 
conventional  optical  coupling  schemes  is  that  the  heat  leak  down  the 
glass  is  negligible.  Based  on  typical  thermal  conductivity  values  of 
glass  (Childs,  at  al . ,  1973),  the  heat  leak  is  estimated  to  be  -  2  ^W. 

In  contrast,  transverse  optical  access  dewars  must  necessarily  allow  a 
path  for  at  least  some  blackbody  radiation  (47  mW/cm  at  300K)  to  enter 
the  sample  area,  as  well  as  the  desired  optical  radiation. 

The  optical  system  is  diagramed  schematically  in  Fig.  6.4.  The 
laser  output  power  could  be  regulated  from  50  mW  to  500  mW.  For  lower 
powers,  a  10  dB  attenuater  was  inserted  into  the  beam  path.  A  glass 
beamsplitter  reflected  about  8  percent  of  the  light  onto  a  silicon 
photodiode  detector.  The  rest  of  the  beam  was  chopped  at  337  Hz  (a 
particularly  low  noise  frequency  for  lock-in  detection),  then  focused 
down  onto  the  core  of  the  optical  fiber.  The  chopper  generated  an 
electrical  signal  which  was  used  to  synchronize  the  two  lock-ins. 
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VII.  THERMOELECTRIC  EFFECT  RESULTS  AND  DISCUSSION 


A.  Voltage  Dependence  Measurements 


A  test  for  the  existence  of  a  thermoelectric  effect  In  tunnel 
junctions  was  first  performed  by  using  the  voltage  measurement  scheme 
described  In  Chapter  II.  The  data  used  for  the  unfolding  analysis  could 
be  used  to  search  for  a  tunneling  thermoelectric  effect.  Illuminating 
the  tunnel  Junction  from  the  lead-bismuth  side  causes  two  observable 
effects  in  the  tunnel  junctions.  The  first  is  a  general  increase  in  the 
tunneling  conductance,  caused  by  the  larger  number  of  excited 
quaslpartlcles  (see  earlier  chapters).  For  constant  current  biasing, 
this  leads  to  a  contribution  to  the  observed  I  dV  |  which  Is  an  even 
function  of  voltage,  and  Is  defined  as  the  symmetric  contribution.  (The 
sign  of  dV  reverses  near  the  origin.)  Of  more  concern  here,  is  the 
thermoelectric  contribution  to  I  dV  |  .  The  thermoelectric  current 
through  the  tunneling  barrier  should  be  only  weakly  voltage  dependent,  as 
Is  clear  from  the  argument  leading  to  Eq .  (5.5).  For  a  fixed  current 
bias,  the  thermoelectric  current,  I^,  must  be  countered  by  a  decrease  In 
the  conventional  quaslpartlcle  current.  This  gives  rise  to  an  observable 
shift  In  the  junction  voltage  given  by: 

dV  -  dV/dl  Ig  (7.1) 


I  . 
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where  dV/dl  is  the  dynamic  resistance  of  the  tunnel  junction. 

Typical  behavior  of  I  dV  |  is  shown  In  Fig.  7.1.  The  absolute 
value  of  dV  Is  plotted  to  show  the  asymmetry,  expeclally  at  low  voltages. 
The  observed  voltage  shift  I  dV(V)  |  could  be  decomposed  into  a  symmetric 
dV^(V),  and  an  antl^syometrlc  dV^(V).  The  values  of  dV^  were  used  in  the 
unfolding  analysis,  and  will  not  be  further  discussed  here.  The  meastired 
anti- symmetric  contribution  Is  plotted  in  Fig.  7.2.  The  falloff  In  dV^ 
for  large  voltages  corresponds  to  the  decrease  in  the  dynamic  resistance 
of  the  tunnel  junction  at  high  voltages,  and  is  expected  from  Eq.  (7.1). 
Fig.  7.3  shows  the  asymmetry  plotted  as  a  function  of  the  dynamic 
resistance  of  the  junction  for  one  tunnel  junction  at  a  fixed  bath 
temperature.  The  straight-*  line  fit,  corresponding  to  a 
voltage-independent  thermoelectric  current  of  *  0.91  nA,  is  good  to 
within  the  accuracy  of  the  data. 

An  estimate  for  Cj^  can  be  obtained  by  comparing  the  observed 
thermoelectric  current  to  theory  (Eq.  (A. 12)).  A  value  of  was 
estimated  by  measuring  the  shift  In  the  Apggj^  ^A1  current  jump 
location  on  the  I(V)  curve.  (See  Chapter  IV.  for  details  on  this  type 
of  calculation.)  The  Inferred  value  of  ^Apgg^  was  then  used  to 
calculate  the  extra  quaslpartlcle  density  in  the  lead-bismuth  according 
to 


FIGURE  7.1.  Observed  laser-induced  voltage  shift  I  dV  |  vs.  V  for 
an  Al-PbBi  tunnel  Junction  illuminated  on  the  PbBi  side.  The 
temperature  was  1.588  K. 


FIGURE  7.3.  Observed  dV  vs.  dV/dl,  the  (voltage  dependent) 
dynaalc  resistance  of  the  tunnel  Junction.  The  straight  line  Is  the 
best  fit  to  a  constant  thermoelectric  current  model. 
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A/A(0)  ■  1  -  2n  (7.2) 

The  altuttinua  film  has  a  relatively  soall  increase  In  the  density  of 
quaslpartlcles  with  sufficiently  high  energy  to  tunnel  (E  ^  For 

the  observed  themoelectrlc  currents  and  the  value  of  Implied  by 

SApb5i»  values  of  0.09  -  0.40  ev“^  were  obtained  for  c^.  These 
estimates  are  an  order  of  magnitude  less  than  the  3.1  eV  ^  predicted  by 
the  one-band,  model  (Eq.  (5.6)),  but  the  direction  of  the  discrepancy  is 
In  qualitative  agreement  with  the  two-band  corrections  of  Gundlach 
(1973). 


B.  SQUIB  Measurements 


The  voltage  measurenent  techniques  are  not  convenient  for  detailed 
measurements  of  the  thermoelectric  effect.  To  measure  the  power  and 
temperature  dependence  of  the  thermoelectric  effect,  the  SQUID  system 
described  la  Chapter  VI.  was  used  to  provide  direct  measurenent  of  the 
zero-voltage  thermoelectric  current,  I^. 

Results  for  a  sweep  In  temperature,  with  constant  (small)  laser 
power,  are  shown  In  Fig.  7.4.  The  sample  temperature  was  Initially  at  3 
K  and  the  temperature  was  reduced  eventually  to  roughly  1.6  K.  The 
current  became  larger  as  T  (and  hence  decreased,  as  expected  from 
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Eq.  (5.16).  There  was  no  noticeable  change  observed  upon  passing 
through  of  the  dirty  aluminum  film.  The  solid  curve  is  proportional 
to  a^(T)~^«  as  my  model  predicts.  The  fit  is  good  for  temperatures  above 
the  lambda  point  of  helltmi.  At  the  lambda  point  there  was  an  abrupt 
decrease  in  the  measured  signal.  This  is  attributed  to  the  improved 
coupling  of  phonons  to  the  now  superfluid  bath,  and  the  resultant 
decrease  in  ‘^eff*  ^  temperature  near  1.85  K  the  signal  disappeared 
to  the  limits  of  resolution  (1  pA)  of  the  experiment.  This  temperature 
corresponded  to  that  at  which  the  Josephson  supercurrent  became  large 
enough  to  be  observed  in  a  separate  experiment  in  which  current  and 
voltage  leads  were  attached  to  the  sample.  Accordingly,  I  interpret  the 
disappearance  of  the  external  thermoelectric  current  as  due  to  its  being 
effectively  shorted  out  by  an  internal  backflow  of  supercurrent.  This 
interpretation  is  supported  by  the  fact  that  at  nonzero  voltages,  where 
there  is  no  dc  supercurrent,  a  thermoelectric  current  ^  measured  in  the 
experiments  reported  in  section  A.,  even  at  temperatures  as  low  as  1.5  K. 

The  SQUID  galvanometer  experiment  was  also  performed  at  fixed 
temperature  while  varying  the  laser  power  P.  Fig.  7.5  shows  the  power 
dependence  of  for  several  temperatures.  As  expected,  the  observed 
currents  become  nonlinear  in  power  for  the  lower  temperatures. 

Fig.  7.6  shows  the  predictions  for  my  simple  model  (Eq.  (5.15)). 
The  value  of  c^^  and  the  value  of  the  phonon  trapping  factor  have  been 
adjusted  to  give  a  best  fit.  For  an  assumed  film  optical  absorptivity  of 
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2  percent,  and  of  4.3  X  10  sec.,  the  fit  value  of  ^  is  15.  As 

before,  the  fit  values  of  c^  (  >0.5  eV~^)  are  an  order  of  magnitude  less 
than  predicted  by  the  one  band  model  calculation  (Eq.  (5.5)),  but  are  in 
qualitative  agreement  with  the  two  band  corrections  proposed  by  Gundlach 
(1973).  While  the  general  magnitude,  shape,  and  sign  of  the  Iq(T)  curves 
are  in  good  agreement  with  the  simple  Rothwarf~Taylor  calculations  done 
here,  the  exact  temperature  dependence  was  less  satisfactory.  No  attempt 
has  been  made  to  include  liquid  helium  effects  which  would  cause  a 
temperature  dependence  to  the  phonon  trapping  factor.  Available 
literature  suggests  that  heat  transport  through  helium  is  strongly 
dependent  on  the  exact  sample  geometry  and  surface  structxure,  and  is 
therefore  difficult  to  calculate  realistically.  The  importance  of  the 
liquid  hellxim  in  determining  phonon  trapping  factor  is  evidenced  by  the 
dramatic  decrease  in  as  the  helium  becomes  superfluid. 

For  a  given  temperature,  the  adequacy  of  the  Rothwarf-Taylor  model 

to  describe  the  energy  relaxation  within  the  superconductor  is  quite 

good.  Fig.  7.7a  shows  the  measured  behavior  of  Iq(P)  in  the-  low 

1/2 

temperature  regime.  The  accuracy  of  the  P  power  dependence  is 

2 

emphasized  in  Fig.  7.7b,  where  Is  plotted  for  the  same  data.  The 

straight-line  behavior  Indicates  that  the  self-consistent  Rothwarf-Taylor 
solution  for  the  number  of  excited  qtiaslpartlcles  fits  the  experimental 
situation  very  well. 
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C.  Comparisons  with  Point  Contact  Experiments 


The  existence  of  a  thermoelectric  current  across  tunneling  junctions 
explains  several  key  dlscrepencles  between  earlier  experiments  and 
existing  theory.  In  particular.  It  Is  useful  to  examine  the  experiments 
of  Clarke  and  Freake  (1972)  In  light  of  temperature  gradients  across 
oxide  barriers.  Their  experiment  consisted  of  a  point  contact  formed  by 
bringing  a  sharpened  Pb  wire  Into  contact  with  a  Pb  foil.  The 
temperatures  of  the  wire  and  the  foil  could  be  Independently  varied  and 
measured.  Although  the  point  contact  supercurrent  effectively  shorted 
out  any  thermoelectric  current.  It  was  possible  to  measure  the  difference 
In  the  magnitude  of  the  critical  current  for  the  two  directions  of 
current  flow.  This  asymmetry  was  Interpreted  In  terms  of  a  thermoelectlc 
quaslpartlcle  current.  Clarke  and  Freake  measured  a  sizeable  current 
whenever  the  point  and  foil  were  at  different  temperatures,  although  the 
size  of  the  effect  varied  over  a  factor  of  50  for  different  point 
contacts.  The  existence  of  any  asymmetry  current  was  surprising  to  those 
who  expected  no  steady-state  thermoelectric  current  to  be  possible  within 
a  superconducting  channel. 

In  an  effort  to  clarify  matters,  Welker  and  Bedard  (1977)  performed 
a  somewhat  similar  experiment  using  Nb  wire.  They  were  extremely  careful 
about  sample  preparation  conditions.  For  very  clean,  oxide-free  point 
contacts,  where  a  small  but  continuous  superconducting  channel  could  be 


assumed  to  exist,  they  found  no  evidence  for  a  thermoelectric  current. 
If  they  allowed  an  oxide  to  form  In  the  point  contact  area,  however,  they 
did  note  a  significant  thermoelectric  current.  The  deduced  current  was 
In  the  same  direction  relative  to  the  temperature  difference  as  was 
observed  by  Clarke  and  Freake,  although  Nb  and  Pb  have  thermopowers  which 
differ  In  sign.  Thus,  the  thermoelectric  voltage  was  of  the  wrong  sign 
to  be  predicted  by  the  thermopower  arguments  of  Clarke  and  Freake,  but  is 
correctly  predicted  by  our  model. 

A  follow  up  experiment  by  Matslnger  ejc  al.  (1978)  was  performed 

using  both  Nb  and  Pb  point  contacts.  Asymmetry  was  looked  for  over  the 

entire  I(V)  curve  of  the  point  contacts,  for  both  directions  of  applied 

temperature  gradient.  For  situations  where  one  electrode  was  raised 

above  Its  T^,  while  the  other  electrode  was  superconducting,  they 

measured  the  expected  normal  thermopower.  With  both  metals  below  T  no 

c 

asymmetry  was  observed  to  the  limit  of  their  sensitivity.  Matslnger,  et 
al .  concluded  that  the  Clarke  and  Freake  measurements  were  "caused  by  a 
surface  layer,”  and  were  not  an  intrinsic  superconducting  effect. 

It  seems  likely  that  the  asymmetric  critical  currents  measured  both 
by  Clarke  and  Freake  and  by  Welker  and  Bedard  are  the  result  of 
thermoelectric  currents  across  oxide  barriers.  The  sign  predicted  here 
for  tunneling  barrier  thermoelectric  currents  agrees  with  experimental 
results.  The  magnitude  of  the  predicted  effect  can  be  calculated  using 
the  point  contact  parameters  given  In  the  Clarke  and  Freake  paper. 
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Integrating  Eq.  (A. 16)  numerically  for  Pb  electrode  temperatures  of  7.0 
and  5.3  K  gives: 

-  (2/eR)  Cj^  (0.08  A^>  (7*3) 

where  A  Is  the  lead  gap  at  5.3  K.  Estimating  Cj^  Is  not  as  simple  as  for 
our  experimental  geometry.  For  Pb  tunnel  junctions,  Basavalah  al . 

(1974)  have  measured  the  oxide  barrier  height  as  1.05  eV.  The  actual 
junction  area  Is  somewhat  poorly  defined  for  a  point  contact.  If  a  1 
jim  area  Is  assumed,  and  the  junction  resistance  Is  1  ohm,  then 
Cj^s  4  eV  ^  and 

^theory  .  o.7  yk  (7.4) 

If  the  junction  effective  area  Is  larger,  the  prediction  Is 

logarithmically  larger.  If  not  all  the  temperature  difference  occurs  at 

the  oxide,  a  smaller  I^  should  result.  Given  the  variance  In  the 

o 

experimental  results,  uncertainties  In  several  of  the  barrier  parameters, 
expeclally  A  and  0,  and  the  tendency  of  the  one- band  approximation  to 
overestimate  c^,  this  Is  In  good  agreement  with  the  plotted  Clarke  and 
Freake  result 

^experiment  .  o.25  yA  (7.5) 


/  , 
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(Note  that  their  experimental  result  has  been  adjusted  by  4  factor  of  two 
because  the  published  values  are  peak-to-peak. )  Coincidentally,  the 
expected  normal  state  thermoelectric  current  for  the  point  contact  la  the 
same  order  of  magnitude,  0.4  ^A.  The  two  Mchanlsms  are  entirely 
Independent;  the  mechanism  reported  here  Is  an  oxide  effect,  dependent 
on  the  barrier  height  and  thickness,  while  the  Seebeck  effect  depends  on 
the  metallic  density  of  states  of  the  electrode  materials. 


V  . 


VIII.  CONCLUSION 


In  suoiaary,  the  study  of  optically  illuminated  Al-PbBl  tunnel 
junctions  has  proven  to  be  effective  in  the  study  of  properties  of 
nonequilibrium  superconducting  films.  A  general  procedure  for  unfolding 
I(V)  measurements  of  asymmetric  tunnel  junctions  has  been  developed  as  a 
tool  for  measuring  the  quasiparticle  energy  distribution.  Absolute 
sensitivity  of  10~^  and  relative  accuracy  of  a  few  percent  are  readily 
obtainable  for  small  perturbations  from  equilibrium.  Results  of 
unfolding  experimental  data  from  optically  illuminated  Al-PbBi  tunnel 

junctions  show  that  the  aluminum  nonequilibrium  state  is  well  described 

* 

by  a  T  model  behavior  over  several  orders  of  magnitude  in  f(E). 

The  enormous  temperature  gradients  present  across  illuminated  tunnel 
junctions  (up  to  10^  K/cm)  lead  to  the  study  of  thermoelectric  effects  in 
tunnel  junctions.  For  nonequilibrium  situations  in  which  a  temperature 
difference  exists  across  a  tunneling  barrier,  I  have  shown  the  existence 
of  a  thermoelectric  current,  which  exists  even  with  no  voltage  difference 
across  the  oxide.  This  current  results  from  the  energy  dependence  of  the 
electronic  tunneling  probability  for  barrier  penetration,  and  is  largely 
independent  of  material  parameters  of  the  electrode  materials.  I  have 
measured  this  thermoelectric  effect  in  an  open  circuit  experiment  as  well 
as  with  a  current  measurement  scheme.  In  both  cases  the  data  are  in  at 
least  qualitative  agreement  with  theoretical  predictions. 


The  existence  of  Che  tunneling  Chenaoelectrlc  current  can  explain 
Che  supercurrent  asynunetrles  studied  in  point  contact  experiments  by  a 
number  of  authors*  The  theory  presented  here  correctly  predicts  the  sign 
and  magnitude  for  these  asymmetries,  as  well  as  the  absence  of  any 
asymmetry  In  very  clean  point  contacts. 

In  addition,  the  temperature  dependence  and  power  dependence  of  the 
thermoelectric  effect  can  be  used  as  a  probe  of  the  nonequilibrium 
superconducting  state.  The  thermoelectric  current  effectively  measures 
the  density  of  excited  quasiparticles.  The  thermoelectric  effect  was 
studied  in  both  the  linear,  weak  perturbation  regime  and  In  the  highly 
non-linear  regime  far  from  equlllbrltim.  In  both  cases  the  behavior  of 
the  superconductor  was  well  described  by  a  self-consistent  solution  to 
the  Rothwarf-Taylor  equations. 

In  conclusion,  the  study  of  the  I(V)  characteristics  of  tunnel 
Juncitons  has  helped  to  provide  insight  into  a  variety  of  nonequilibrium 
superconducting  phenomena^  It  is  my  hope  that  the  tools  developed  here 
will  find  future  uses  in  exploring  the  superconducting  state. 


o  o  o  o 
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APPENDIX  I.  UNFOLDING  PROGRAM 


C 

C 

C 

C 

C 

C 

C 

C 

C 
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35 

36 
34 

52 


51 

50 


3 

1 

6 


C 
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ASYMMETRIC  JUNCTION  UNFOLDING  PROGRAM 

WRITTEN  BY  ANDY  SMITH 

LINES  BEGINNING  WITH  C  ARE  CGHMENTS 

LINES  BEGINNING  WITH  D  ARE  IGNORED  EXCEPT  WHILE  DEBUGGING 
COMMDN  /DEBUG/IDBG(2) 

COMMON  Y(IOO) 

DIMENSION  C0EFFS(20),VEES(100),C(21) 

DIMENSION  0BS(21, 100) ,A(20) ,STDA(20) 

INPUT  FORMAT 
LINE  NO. 

1  BLANKS  UNLESS  DEBUGGING 

2  LARGER  GAP,  SMALLER  GAP  (MEV),  TEMPERATURE  (K), 
RESISTANCE  (OHMS)  ALL  FLOATING  POINT  NUMBERS 

3-8  ENERGY  BIN  LIMITS  (21  VALUES,  4  PER  LINE) 

9-  VOLTAGE  (MV), CURRENT  (MA)  DATA,  1  PAIR  PER  LINE 

OUTPUT  IN  FILE  1  CONTAINS  BEST  FIT  F-TWIDDLE  VALUES, 
UNCERTAINTIES  IN  FIT  VALUES,  AND  FINALLY  DATA  WITH  FIT 
COMMON  /PRAM/V, TEMP, DELTAA, DELTAS, E(21) 

READ(20,110)IDBG 

F0RMAT(5I1) 

READ ( 20 , 2)DELTAS ,DELTAA,TEMP , RESIST 

WRITE(1,2)  deltas, DELTAA, TEMP, RESIST 

R£AD(20,34)E 

WRITE(1,35)E 

FORMAT(4F7.3) 

DO  36  K-1,21 
E(K)»E(K) /DELTAS 
F0BMAT(4F) 

DO  1  NOBS-1,100 
READ(20, 2 ,END-6 )V ,CURENT 
CURENT-CURENT*RESI ST 
IF(V.LT.DELTAS-DELTAA)G0  TO  50 
TYPE  51, NOBS 

FORMAT('  ATTEMPT  TO  USE  VOLTAGE  TOO  HIGH  ON  OBS  ',15) 

GO  TO  52 
VEES(NOBS)^ 

ENTER  DATA  AND  G  VALUES 
CALL  CRAMK2(C) 

DO  3  KVAR-1,20 
OBS(KVAR,NOBS)-C(KVAR) 

0BS(21,N0BS)-CURENT 

CONTINUE 

NOBS-NOBS-1 

WRITE(22, 37) ( ((OBS(MM,M)  ,MW1 , 21)  ,VEES(M)  )  ,M-1  ,NOBS) 

IS  THERE  ENOUGH  DATA? 
tF(N0BS.GE.21)G0  TO  55 
TYPE  56,NQBS 

PORMAT('  ONLY  ',15,'  DATA  POINTS  ') 


no  o  o  r>  o 


102 


STOP 

55  CONTINUE 

2  F0RMAT(4F) 

MAYBE  I  SHOULD  CHECK  FOR  FIT  WITH 
A  PARALLEL  CONDUCTANCE  AS  ONE  OF  THE  BINS 
IF  SO,  REPLACE  LAST  BIN  BY  CONDUCTANCE 
IF(  IDBG(2)  .NE.  1)  GO  TO  120 
DO  121  K-1,N0BS 
121  0BS(20,K)-VEES<K) 

120  CONTINUE 
37  F0RMAT(6E13.4) 

DO  THE  JOB 

CALL  FITIT(0BS,N0BS,21,Y,C0EFFS,STDA) 

WRITE( 1 , 10)C0EFFS , STDA 
CHI-0. 

DO  77  K-1,N0BS 

CHI-CHl+( OB  S ( 2 1 , K)-Y ( K) ) **2 /OBS ( 2 1 , K) **2 
77  WRITE(1,10)VEES(K),0BS(21,K),Y(K) 

STDIffiV-SQRT(CHI/FL0AT(N0BS-21) ) 

WRITE (1,10)  STDDEV 
10  FORMAT (2X,5E1 1.4) 

25  F0BMAT(1X,5E12.3) 

27  F0RMAT(5E15.6) 

STOP 

END 

SUBROUTINE  CRANK2(OUr) 

AT  A  GIVEN  VOLTAGE,  CALL  ROUTINES  TO  CALCULATE  ALL  OF  THE 
G  VALUES 

COMMON  /PRAM/VREAL , TEMP , DELTAA, DELTAS , E ( 2 1 ) 

COMMON  VD,BD,DADB2,X,Y 
DIMENSION  OUT (11  ) 

BD-1 1. 6*DELTAS/TEMP 
VD-VREAL/DELTAS 
DADB2-(DELTAA/DELTAS) **2 

our(2i)-o 

999  F0FMAT(I4) 

DO  3  K-1,20 
8811  F0RMAT(5F10.3) 

OUT(K)-DELTAS*CCC(E(K) ,E(K+1 ) ) 

3  OUT(21)-ODT(21)+OUT(K) 

CONTINUE 

2  CONTINUE 

4  F0RMAT(11E10.3) 

RETURN 

STOP 

END 

FUNCTION  CCC(EMIN,EMAX) 

C  COMPUTE  G  VALUE  FOR  A  GIVEN  ENERGY  RANGE 
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EXTERNAL  CO 

C  CONTRIBUTION  FROM  QUASIPARTICLES  WITHIN  ENERGY  RANGE 

COMMON  VD,BD,DADB2,KI,T2ERO 
DIMENSION  P0INT(4) 

101  F0RMAT(F10.3,I5,F10.3) 

TOTAL -0 . 

ITYPE-0 

100  FORMAT(5F10.3,I5,F10.3) 

888  F0RMAT(I5) 

C  TEST  LIMITS  OF  INTEGRATION 

P0INT(1)-1.-VD 
P0INT(2)-1. 

P0INT(3)-1.+VD 

P0INT(4)-AMAX1(1 .+8. /BD ,P0INT(3)+1 . ) 

DO  1  K-1,3 

C  DIVira:  INTO  REGIONS 

Kl-K 

El-AMAXl (POINT(K) ,EMIN) 

103  FORMAT(5F10.4) 

C  AVOID  SINGUURITIES 

1F(E1.GT.EMAX)G0  TO  99 
IF(E1.GT.P0INT(K+1))G0  TO  1 
E2-AMIN1 (EMAX , P0INT(K+1 ) ) 

T1-SQRT(E1-P01NT(K) ) 

T2-SQRT(E2-POINT(K) ) 

TZERD-POINT(K) 

ADDIN-SIMPSN(C0,T1 ,T2 ,22) 

TOTAL-TOTAL+ADDIN 

IFCITYPE.EQ.DTYPE  106,  VD,KI, El  ,E2  ,T1  ,T2 ,  ADDIN 
106  F0RMAT(F5.2,I5,4F6.3,E10.3) 

1  CONTINUE 

99  CCC-2.*T0TAL 

RETURN 
END 

FUNCTION  C0(T) 

C  0ADB2<1 . 

COMMON  VD,BD,DADB2,K1,TZER0 
C  ALL  ENERGIES  EXPRESSED  FOR  LARGER  DELTA 

ENT(X)-X/SQRT(X+1 . ) 

EN2 (X)-X/S(SIT<X**2-DADB2) 

ENl (X)-X/SqRTCX**2-l . ) 

E-T**2-KrZER0 

102  F0RMAT(3F10.4,I6) 

CO-0. 

C  SWITCH 

GO  TO  (1,2,3),KI 
3  CO— ENT(E“VD)*EN2(E)+ 

+  •r*ENl(E)*(EN2(E4VD)-EN2(E-VD>)+ 


o  o  o  o  o  o  o 
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+  T*EN1(E+VD)*EN2(E) 

GO  TO  9 

2  CO-ENT(E)*(EN2(E-WD)-EN2(E-VD))+ 

+  T*EN1(E4VD)*EN2(E) 

GO  TO  9 

1  C0-C0+-ENT(E+VD)*EN2(E) 

9  CO-CO/ (l.+EXP(E*BD)) 

RETURN 
END 

THIS  PROGRAM  IS  CALLED  BY  THE  MAIN  PROGRAM 
GIVEN  DATA,  IT  LOADS  THE  DATA  AND  CROSS-CORRELATES 
IT  THEN  CALLS  A  MATRIX  INVERSION  SUBROUTINE  TO  UNFOLD 
SUBROUTINE  FITIT(OBS,NOBS,NVAR,Y, A.STDA) 

C»!MON  /DEBUG/IDEBUG 
DOUBLE  PRECISION  ALPHA,DET,WA, WB 
DIMENSION  ALPHA(20, 20) ,BETA(20) ,ST0RE(20, 20) 

DIMENSION  SH0WIT(20),0UTPT(20) 

DIMENSION  WA(35),WB(35) 

DIMENSION  0BS(NVAR,N0BS) ,A(20) , STDA(20) ,Y(N0BS) 
WGT(I0BSVR)-1 . /OBS(NVAR, I0BSVR)**2 
MVAR-NVAR-1 
TYPE  60 

60  FORMATC'  ENTERRED  FITIT  ') 

DO  1  KV-1,MVAR 
DO  2  LV-KV,MVAR 
SUM-0 

CALCULATE  CROSS  CORRELATION  MATRIX 
DO  3  IIOBS-l.NOBS 
lOBS-IIOBS 

3  SOM-SUM+  OBS(KV,IOBS)*OBS{LV,IOBS)*WGT(IOBS) 
ALPHA(KV,LV)-SUM 

2  ALPHA(LV,KV)-SUM 
1  CONTINUE 

WRITE(23,55)  OBS 
WRITE(23,55)  ZERO 
WRITE(23,55)ALPHA 
TYPE  61 

C  SET  UP  BETA  COEFFICIENTS 

61  FORMATC'  SETTING  DP  BETA  COEFFICIENTS') 

DO  4  KV-1,MVAR 

SUM-O. 

DO  5  IIOBS-l.NOBS 
lOBS-IIOBS 

5  SUM-SUM-K)BS (KV , lOBS ) *OBS(NVAR, lOBS ) *WGT (lOB  S ) 

4  BETA(KV)-SUM 
DO  6  KA-1,MVAR 
DO  6  KB-1,MVAR 

6  STORE(KA,KB)-ALFHA(KA,KB) 


o  o  o  o  t» 
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C  CALCULATE  EPSILON  MATRIX 

TYPE  63 

63  FORMAT('  INVERTING') 

CALL  MATINV (ALPHA, MVAR,DET,WA,WB) 

D  WRITE(23,55)ALPHA 

ZERO-0. 

WRITE(23,55)  ZERO 
WRITE(23,55)  BETA 
55  FORMAT ( IX, 5E1 1.3) 

DO  8  KVAR-1,MVAR 
A(KVAR)-0 
DO  8  LVAR-l.MVAR 

8  A(KVAR)-A(KVAR)+ALPHA( KVAR, LVAR) *BETA(LVAR) 

SUMB-0. 

DO  10  I0BS-1,N0BS 
SUM-  .0 

DO  9  IVAR-1,MVAR 

9  SUM-SUM+A(IVAR)*OBS(IVAR,IOBS) 

Y(IOBS)-SUM 

10  SUMB-SUMB+(SUM-OBS(NVAR, lOBS) )**2*WGT(IOBS) 

VARNCE-SUMB/ FLQAT(NOB S-NVAR) 

THIS  IS  THE  NEW  STANDARD  DEVIATION  SECTION 

DO  200  KVAR-1,MVAR 
STDA(KVAR)-0. 

DO  202  I0BS-1,N0BS 
SUM-0. 

DO  201  LVAR-l.MVAR 

201  SUM-SDM+ALPHA(KVAR, LVAR) *OBS (LVAR, lOBS) 
SUM-SUM*WGT(IOBS) 

C  SUM  IS  NOW  DU,  THE  MATRIX  RELATING  F  TO  I 

STDA(KVAR)-STDA(KVAR)+(SDM* (Y (lOBS )-OBS(NVAR, lOBS ) ) )**2 

202  CONTINUE 

200  STDA(KVAR)-SQRT(STDA(KVAR)) 

70  CONTINUE 

TYPE  lll.IIEBDG 

111  F0RMAT('  lUEBUG  IS  ',15) 

IF  (IDEBUG  .NE.  1)RETURN 
C  NORMALIZE  CROSS-CORRELATIONS 

WRITE(24,100) 

100  F0RMAT('  NORMALIZED  CROSS  CORRELATION  ') 

101  FORMAT(X,21F6.2) 

DO  102  K-1,20 
DO  103  L-1,20 

103  OUTPT(L)-STORE(K,L)/SQRT(ABS(STORE(K,K)*STORE(L,L))) 

102  WRITE(24,101)  OUTPT 
WRITE(24,123) 


123  FORMAT('  FRACTIONA.L  CONTRIBUTIONS  ') 

C  FRACTIONAL  COHRIBUTIONS 

DO  104  llOBS-l.NOBS 
lOBS-IIOBS 
DO  105  K«l,20 
SDM-0. 

DO  106  L-1,20 

106  SUM-SUM+ALPHA(K,L)*WGT(10BS)*0BS(L,I0BS) 
105  0UTPT(K)-SUM/A(K)*0BS(21,I0BS) 

104  WSITE(24,101)0UTPT 

72  F0KMAT(15,10F7.4) 

RETURN 

END 

SUBROUTINE  MAT1NV{A,N,D,L,M) 

C  MATRIX  INVERTER 

C  ARGUMENTS  IN  ORDER  ARE: 

C  INPUT  MATRIX  (DESTROYED  IN  COMPUTATION) 

C  DIMENSIONALITY  OF  MATRIX 

C  DETERMINANT  OUTPUT 

C  WORK  VECTOR  LENGTH  N 

C  WORK  VECTOR  LENGTH  N 

C  MODIFIED  FOR  DOUBLE  PRECISION 

IMPLICIT  DOUBLE  PRECISION  (A-«,0-Z) 
DIMENSION  A(1),L(1),M(1) 

D-1. 

NK— N 

DO  80  K«l,N 

NK-NK+N 

L(K)-K 

M(K)-K 

KK-NK+K 

BIGA-A(KK) 

DO  20  J-K,N 
IZ-N*(J-1) 

DO  20  I-K,N 
IJ-IZ+I 

10  IF<DABS(BIGA)-DABS(A(IJ)))  15,20,20 
15  BIGA-A(IJ) 

L(K)-I 

M(K)-J 

20  CONTINUE 
J-L(K) 

IF(J-K)  35,35,25 
25  KI-K-M 

DO  30  I-1,N 
Kl-KI+N 
HOLD— A(KI) 

JI-KI-K+J 
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A(KI)-A(JI) 

30  A(JI)-H0LD 
35  I-M(K) 

IF(I-K)  45,45,38 
38  JP-N*(I-1) 

DO  40  J-1,N 
JK-NK+J 
JI^P+J 
HOLD— A(JK) 

A(JK)-A(JI) 

40  A(JI)-HOLD 

45  IF(BIGA)  48,46,48 

46  D*0. 

RETURN 

48  DO  55  I-1,N 

IF(I-K)  50,55,50 
50  IK-NK+I 

A(IK)-A(IK) /BIGA*-1. 

55  CONTINUE 

DO  65  I»1,N 

IK-NK+I 

HOLD-A(IK) 

IJ-I-N 
DO  65  J-1,N 
IJ-IJ+N 

IF(I-K)  60,65,60 
60  IF(J-K)  62,65,62 

62  KJ-IJ-I-HC 

A(IJ)  -H0LD*A(KJ)  +A(IJ) 

65  CONTINUE 

C  DIVIDE  ROW  BY  PIVOT 

KJ-K-N 
DO  75  J-1,N 
KJ-KJ+N 

IF(J-K)  70,75,70 
70  A(KJ)-A(KJ)/BIGA 

75  CONTINUE 

C  PRODUCT  OF  PIVOTS 

D-D*BIGA 

C  REPLACE  PIVOT  BY  RECIPROCAL 

A(KK  )-l./BIGA 
80  CONTINUE 

C  FINAL  ROW  AND  COLUMN  INTERCHANGE 

K-N 

100  K-(K-l) 

IF(K)150,150,105 
105  I-L(K) 

IF(HC)  120,120,108 


/ 
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JQ-N*(K-1) 

JR-N*(I-1) 

DO  110  J-1,N 

JK-JQ+J 

HOLD-A(JK) 

JI-JR+J 
A(JK)— A(JI) 

A(JI)  -HOLD 
J^(K) 

IF(J-K)  100,100,125 

KI-K-N 

DO  130  I-1,N 

KI-KI+N 

HOLD-A(KI) 

A(KI)  — A(JI) 

A(JI)  -HOLD 
GO  TO  100 
RETURN 
END 

FUNCTION  SIMPSN(MHAT,XMIN,XMAX,N) 

NUMERICAL  INTEGRATION  ROUTINE  USING  SIMPSONS  RULE 
DX-(XMAX-XMIN)*.5/N 
X-XMIN40X 

SUM-WHAT ( XMIN ) - WHAT ( XMAX) 

DO  1  K-1,N 

SUM-SUM+4 . *WHAT(X)+2 . *WHAT(X+DX) 

X-DX*2.+X 

SIMPSN-DX*SUM/3. 

RETURN 

END 
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APPENDIX  II.  THERMOELECTRIC  TUNNELING 

A.  Tunneling  Hamiltonian 

To  calculate  correctly  the  current  through  a  superconducting  tiinnel 
junction,  it  is  necessary  to  take  account  of  the  coherence  factors  which 
describe  the  quasiparticle  wave  functions.  Barrier  transmission 
probability  amplitudes  "Vcq  describe  coupling  between  states  with  momenta 
k  in  one  metal,  and  q  in  the  other  metal.  For  weak  coupling  through  the 
oxide  barrier,  we  can  assume  the  familiar  tunneling  hamiltonlan  of  the 
form 


I 

k,q 


^kq  4  \  *  \k  ‘^q 


(A.l) 


where  c^  annihilates  an  electron  from  the  q  side  of  the  junction.  The 
sum  is  taken  over  electron  spins  as  well  as  momenta.  The  probability 
amplitudes  will  be  seen  to  depend  on  the  type  of  tunneling 

(electron- like  or  hole-like)  in  addition  to  the  initial  and  final  states, 
k  and  q. 


The  tunneling  hamiltonlan  must  be  rewritten  in  terms  of 
quasiparticle  states  in  order  to  compute  the  allowed  tunneling  current  in 
the  superconducting  case.  Following  the  notation  of  Tinkham  (1972) 
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'^ekO'*’  ''k\kl^^“q  '^eqO'^q  hql^ 

■*■  ^qk^“k  "^ekO"**  ^k^hkl^^“q  ^eqO"^  ^q  ^hql)]  ^  ^ 

The  hamiltonlan  consists  o£  eight  terms,  uhlch  fall  Into  two 
categories.  There  are  four  terms  Involving  one  creation  and  one 
annihilation  operator.  These  terms  describe  tunneling  of  quaslpartlcles 
from  one  side  of  the  film  to  the  other.  The  other  four  terms  each  have 
either  two  creation  or  two  annihilation  operators,  and  describe 
pair-breaking  or  pair-forming  tunneling  processes. 

The  effects  of  these  terms  are  listed  In  table  A.l.  In  describing 
the  energies  and  E^,  the  excitation  picture  has  been  used,  so  that  the 
energies  are  all  expressed  as  positive  numbers  relative  to  the  Fermi 
surface.  Note  that  for  a  given  energy  E^,  tunneling  will  not  necessarily 
be  allowed  in  a  given  channel  because  of  the  existence  of  a  gap  In  the 
density  of  states.  (Tunneling  is  disallowed  If  E^  <  Each 

tunneling  term  has  been  classified  as  electron-like  (e) ,  hole-like  (h) , 
or  pair-breaking  (pb).  The  sense  of  the  voltage  across  the  tunnel 
junction  has  been  defined  such  that  electrons  gain  energy  In  tunneling 
from  the  k-slde  film  to  the  q-side  film. 

For  simplicity  we  will  restrict  our  attention  to  the  case  that 
eV  <  +  ^q^*  this  situation  pair-breaking  tunneling  is 


Term  in Probability  E  Electron  Type  of  Mean  Electron 

Factor  **  Transfer  Tunneling  Energy  in  the 
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prohlbltted  by  conservation  of  energy. 

By  Feiml's  Golden  Rule,  the  quasiparticle  current  is  simply  (Smith, 

1975) 

Iqp  -  /dE  i  I  ^  N^(E)Nq(E+eV) [f^^(E)-fq(e+eV)] 

-uj  «q  \(E)Nq(E-eV)(f^(E)-fq(E-eV)I)  (A. 3) 

The  sum  is  over  both  excitation  branches  for  both  films  (four  terms). 

The  sum  over  spins  contributes  a  factor  of  2,  which  has  been  incorporated 

2 

into  the  equation.  The  terms  involving  u.  u“  describe  the  probability 

K  q 

of  tunneling  of  electron  excitations  from  the  k*>side  film  to  the  q-side 

2  2 

film,  shown  as  A-type  processes  in  Fig.  5.1.  The  v.  v  terms  involve 

K  q 

transferring  charge  in  the  other  direction  (i.e.  B^type  hole  tunneling). 

B.  Zero  Voltage  NS  Tunneling 

To  fxurther  simplify  the  calculations,  we  will  study  the  special  case 
of  zero-voltage  current  for  a  normal  metal -Insulator- superconductor  (NS) 
tunnel  junction  (V-0  and  Aq“0)-  Under  these  restrictions,  Eq .  (A. 3) 

may  be  reduced  by  performing  the  sum  over  the  electron  and  hole  branches 
of  the  normal  metal  using  the  values: 
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q  <  <1f 

q  >  qp  (A. 4) 

q  <  q? 
q  >  qp 

Therefore 

m 

-  4n  e/ti£  fdE  (  1  Tj^q  1  -  I  |  ^u^)N^(E)  tf,^(E)-fq(E)l 

(A. 5) 


The  sum  now  Is  only  over  the  lt>kp  and  k<kj.  branches  of  the 

2 

superconductor.  As  was  the  case  before,  the  u^  terms  describe  tunneling 
from  electron-like  states  in  the  superconductor,  through  the  oxide,  and 
into  the  normal  metal  with  an  energy  E  above  the  Fermi  surface.  For  this 
electron  tunneling  it  will  be  assumed  that  to  lowest  order  in  the  energy 
(measured  from  the  Fermi  energy) 

T  2 

-  1  +  C,  E  q  >  q-  (A. 7) 

T(0) 


where  T(0)  is  the  tunneling  probability  amplitude  at  the  Fermi  surface. 

2 

Similarly  for  hole  tunneling  (v^^  terms) 
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» 


T(0)i 


q  <  qp  (A. 8) 


Substituting  these  values  for  ^kq  Into  Eq.  (A. 5)  yields 

« 

-  (4n  eM)  I  T(0)  |  [^dE  E  (vj  +  t4)\(E)  [fj^(E)-fq(E)  ] 

(A. 9) 


2  2 

Using  the  Identity  (v^  +  u^  *  1),  we  can  eliminate  the  dependence  on 
coherence  factors  In  Eq.  (A. 9).  The  siim  over  the  k>kp  and  k<kp  branches 
then  simply  Introduces  a  factor  of  two,  leaving: 


iNS  .  2  [(4we/h)  I  T(0)  I  {  dE  E  N^(E)  [f,^(E)-f  ^^(E)  ]  (A. 10) 

The  factor  [(4iTe/fi)  1  T(0)  |  may  be  replaced  by  (1/eR). 

Several  limiting  cases  of  Eq.  (A. 10)  are  worthy  of  note.  If  both 
electrodes  are  at  the  same  tenperature,  there  will  be  no  net  current 
flow.  On  the  other  hand.  If  the  q-slde  film  Is  at  zero  temperature 
(fq  "0),  Eq.  (A. 10)  reduces  to  the  not  too  surprising  result 


(  , 


i 
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I  -  (2/eR)  (  dE  E  N^(E)  fjj_(E)  (A. 11) 


-  (2/eR)  A  <S.^> 


(A. 12) 


where  Is  the  average  qtiaslparticle  energy  In  the  k-slde  film  and  n^ 

Is  the  normalized  quaslpartlcle  density. 


If  the  quaslpartlcle  occupation  of  the  k-side  film  Is  a  Fermi-Dlrac 
distribution  with  temperature  and  the  q^-slde  film  has  temperature 
then  the  Integral  of  (A. 10)  can  be  expressed  as  a  sum  of  modified  Bessel 
functions. 


-  (2/eR)Cj^A^Z  (-l)“^^fKo(»A/k\)+(kT,^/fflA)Ki(mA/kT^) 

n^l 


-Ko<mA  /kTq)-(kTq/mA  )Ki  (mA  ]  (A.  13) 


The  series  converges  rapidly  at  low  temperatures.  The  leading  low 
temperature  behavior  Is 

Mc  1/2  ^^^k  1/2  A /kX_ 

’  [(kT./A)^'^e  -  (kT  /A)^'  e  (A.14) 
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C.  SIS  Tunneling 

The  more  general  case  of  tunneling  between  two  superconductors  can 
be  calculated  In  a  manner  similar  to  that  used  to  derive  Eq.  (A. 10). 
Care  should  be  taken  to  recognize  the  electron/hole  dependence  of 

A  K 

Denoting  the  two  matrix  elements  as  T^^  and  T^j^  respectively,  Eq. 
(A. 3}  becomes 

Iqp  -  4«e/h£  ^dE  (  I  T^q  I  ^  vJv^N^^(E)Nq(E+eV) [f^^(E)-fq(E+eV) ] 

■'  ^kq*  ^  ‘4uq\(E)Nq{E-eV)Ifj^(E)-fq(E-eV)]}  (A. 15) 

Summing  over  the  branches  yields 

^qp  -  e/h  i  dE  1  1  Tjjq  I  ^  \(E)N^(E+eV)[fj^(E)-fq(E+eV)l 
-I  Tjql  ^  N^(E)Nq(E-eV)[f^(E)-fq(E-eV)]}  (A. 16) 

Using  the  oxide  electron  energies  given  in  the  previous  table,  this  can 
be  approximated  as 

« 

Iqp  -  (1/eR)  /dE{[l+Ci<E+eV/2)]N^(E)Nq(E+eV)[f^^(E)-fq(E+eV)] 


i 


-t  l-c^  (E-eV/2  )  ]N^(E  )Nq(E-eV)  { f j^(E )-f q(E-eV)  1 }  (A .  1 7) 


The  eV/2  terms  account  for  energy  gain  (or  loss)  of  the  electron  as  it 
traverses  the  oxide  barrier  In  the  presence  of  an  electric  field. 
Generally  the  thermoelectric  constrlbutlon  to  the  current  is  only 
slightly  voltage  dependent.  The  coefficient  c^  Is  only  weakly  voltage 
dependent  (Gundlach,  1973)  and  can  generally  be  treated  as  a  constant  for 
low  voltages. 

In  all  cases  examined  here,  the  basic  temperature  dependence  Is  well 
described  by  the  approximation  made  In  Eq.  (5.10). 


/ 


118 

REFERENCES 

A.  G.  Aronov  and  B.  Z.  Splvak,  Sov.  J.  Low  Temp.  Phys  641  (1978). 

O.  F.  Balkashln,  I.  K.  Yanson,  and  A.  V.  Khotekevich,  Flz.  Nizk.  Temp. 

28  (1979)  [Sov.  J.  Low  Temp.  Phys.  12  (1979)]. 

S.  Basavaiah,  J.  M.  Eldrldge,  and  J.  Matlsoo,  J.  Appl.  Phys.  45,  457 

(1974). 

P.  Bevlngton,  Data  Reduction  and  Error  Analysis  for  the  Physical 
Sciences,  (McGraw-Hill,  New  York,  1969). 

J.  J.  Chang  and  D.  J.  Scalapino,  Phys.  Rev.  Lett.  37,  522  (1976). 

J.  J.  Chang,  Phys.  Rev.  Lett.  39,  1352  (1977). 

J.  J.  Chang  and  D.  J.  Scalapino,  J.  Low  Temp.  Phys.  31,  1  (1978). 

J.  Clarke,  B.  R.  FJordbdge,  and  P.  E.  Llndelof,  Phys.  Rev.  Lett.  642 

(1979). 


J.  Clarke  and  S.  M.  Freake,  Phys.  Rev.  Lett.  588  (1972). 


119 


J.  Clarke  and  M.  Tlnkham,  Phys.  Rev.  Lett.  106  (1980). 

R.  C.  Dynes,  V.  Narayanamurtl,  and  J.  P.  Gamo,  Phys.  Rev.  Lett.  1509 
(1978). 

W.  Elsenmenger,  K.  Lassman,  H.  J.  Trumpp,  and  R.  Krauss,  Appl.  Phys.  12, 
163  (1977). 

H.  Ekrut  and  A.  Hato,  J.  Appl.  Phys.  1686  (1980). 

Yu  M.  Gal'perln,  V.  L.  Gurevich,  and  V.  I.  Kozub,  Zh.  Eksp.  Teor.  Fir. 
1387  (1974)  [Sov.  Phys.-JETP  W,  680  (1974)]. 

J.  C.  Garland  and  D.  J.  Van  Harlingen,  Phys.  Lett.  47A,  423  (1974). 

Yu.  M.  Gerbshteln,  E.  I.  Nikulin,  and  B.  Z.  Splvak,  Fir.  Tverd.  Tela  20, 
1067  (1978)  [Sov.  Phys.  Solid  State  616  (1978)]. 

K.  E.  Gray  and  U.  W.  Wlllemsen,  Bull.  Am.  Phys.  Soc.  V^,  328  (1979). 

K.  E.  Gray  and  H.  W.  Willemsen,  J.  Phys.  Colloq.  39,  C6/513  (1979). 


\ 


K.  R.  Gundlaeh,  J.  Appl.  Phys.  5005  (1973). 


120 


W.  A.  Harrison,  Phys.  Rev.  123,  85  (1961). 

T.  E.  Hartaan,  J.  Appl.  Phys.  35,  3283  (1964). 

D.  F.  Heidel  and  J.  C.  Garland,  Bull.  Am.  Phys.  Soc.  2£,  411  (1980). 

P.  Hu,  R.  C.  Dynes,  and  V.  Narayanamurti,  Phys.  Rev.  B  10,  2786  (1974). 

S.  B.  Kaplan,  J.  Low  Temp.  Phys.  37,  343  (1979). 

1 

S.  B.  Kaplan,  D.  D.  Chi,  D.  N.  Langenberg,  J.  J.  Chang,  S.  Jafarey,  and  ^ 

D.  J.  Scalaplno,  Phys.  Rev.  B  14,  4854  (1976). 

S.  B.  Kaplan,  J.  R.  Klrtley,  and  0.  N.  Langenberg,  Phys.  Rev.  Lett.  22.* 

291  (1977). 

J.  R.  Klrtley,  D.  S.  Kent,  S.  B.  Kaplan,  and  0.  N.  Langenberg,  J.  Phys. 

Colloq.  39,  C6/511  (1978). 

M.  K.  Konkin  and  J.  G.  Adler,  J.  i^pl.  Phys.  50,  8125  (1979). 

4 


D.  N.  Langenberg,  Proc.  Low  Temp.  Phys.  LT-14  223,  (North-Holland 

Publ.  Co.,  Amsterdam,  1975). 


4 


D.  N.  Langenberg,  In  Feskorperprobleme  XIV,  67,  H.  J.  Quelsser  (ed.)>  (F. 
Vleweg  and  Sohn,  Braunachvieg ,  Germany,  1974). 

A.  A.  J.  Matsinger,  R.  de  Bruyn  Ouboter,  and  H.  van  Beelen,  Physica  93B, 
63  (1978). 

R.  G.  Melton,  J.  L.  Paterson,  and  S.  B.  Kaplan,  Phys.  Rev.  B  2^,  1858 

(1980). 

C.  S.  Owen  and  D.  J.  Scalaplno,  Phys.  Rev.  Lett.  1559  (1972). 

W.  H.  Parker,  Phys.  Rev.  B  12,  1667  (1975). 

W.  H.  Parker  and  W.  D.  Williams,  Phys.  Rev.  Lett.  W,  924  (1972). 

W.  H.  Parker,  Phys.  Rev.  B  3667  (1975). 

C.  J.  Pethick  and  H.  Smith,  Phys.  Rev.  Lett.  43,  640  (1979). 

A.  Rothwarf  and  B.  N.  Taylor,  Phys.  Rev.  Lett.  27  (1967). 

A.  Schmid  and  G.  Schdn,  Phys.  Rev.  Lett.  973  (1979). 

1.  K.  Schuller  and  C.  M.  Falco,  Fix.  Nizk.  Temp..  939  (1978)  [Sov.  J. 


/ 


122 


Low  Temp.  Phys.  446  (1978)]. 

A.  D.  Smith,  M.  J.  Skocpol,  and  M.  Tinkham,  Phys.  Rev.  B  21,  3879  (1980). 

L.  N.  Smith,  Fh.  D.  thesis,  unpublished  (1975). 

M.  Tinkham,  Phys.  Rev.  B  1747  (1972). 

C.  J.  Van  Harlingen,  Bull.  Am.  Phys.  Soc.  411  (1980). 

D.  J.  Van  Harlingen  and  J.  C.  Garland,  Solid  State  Commun.  419 

(1978). 

N.  K.  Welker  and  F.  D.  Bedard,  in  SQUID:  Superconducting  Quantum 
Interference  Devices  and  Their  Applications,  H.  D.  Hahlbohm  and  H.  Liibbig 
(eds.),  (Walter  de  Gruyter,  New  York,  1977),  200. 

H.  W.  Willemsen  and  K.  E.  Gray,  Phys.  Rev.  Lett.  41,  812  (1978). 


H.  W.  Willemsen  and  K.  E.  Gray,  Bull.  Am.  Phys.  Soc.  24,  329  (1979). 


123 


ACKNOVILEDGEMENTS 

It  is  a  pleasure  to  thank  Profs.  M.  Tlnkham  and  W.  Skocpol  for 
their  support  and  encouragement  throughout  -  the  course  of  these 
experiments.  Their  criticisms  and  suggestions  during  manuscript 
preparation  have  been  both  incisive  and  numerous.  Discussions  with  Drs. 
T.  Klapwijk  and  L.  Smith  were  invaluable  in  separating  wheat  from  chaff 
in  theoretical  models  to  explain  observations. 

I  am  grateful  to  the  Gordon  McKay  graduate  students  and  staff 
without  whose  technical  and  moral  support  this  project  would  have  been 
impossible.  G.  Blonder  and  D.  Frank  deserve  credit  for  construction  of 
Che  data  acquisition  computer.  I  would  also  like  to  thank  L.  DeFeo, 
whose  ability  to  build  the  experimental  equipment  far  exceeded  my  ability 
to  draw  it. 

This  work  was  supported  in  part  by  the  National  Science  Foundation 
and  the  Joint  Services  Electronics  Program. 


DlStHiauTlON  LIST  FOR  ONR  ELECTROMIC  KHO  SOLID  STME  SCIEWES 


Director 

Advanced  Raa*arch  Rrojacta  Aigancy 
Attnt  Tachnical  Library 
1400  Hilaon  Boulavard 
Arlington,  Virginia  23200 


OfCicd  of  Naval  Reaearch 
electronics  Prograa  Office  (Code  427) 
800  North  (Quincy  Street 
Arlington,  Virginia  22217 


Office  of  Naval  Research 
code  lOS 

800  North  Quincy  Street 
Arlington,  Virginia  22217 


Director 

Naval  Reaearch  Laboratory 
45S5  Overlook  Avenue.  S.w, 
Nashingtofl.  D.C.  2037S 
Attn:  ftchnical  Library  (S  c(>a) 
Code  $200  (1  copy) 

$210  <1  copy) 

5270  <i  copy) 

8400  (1  copy) 


Office  of  the  Director  of  Defense 
Research  and  Engineering 
Office  of  the  Assistant  Dlr<*ctcr 
Electronics  4  Physical  Sciences 
The  Pentagon,  Room  301079 
ttashington,  DC  20301 


Defense  Docuawntation  Center  <12  cps) 
Ceaeron  Station 
Alcsandrla.  Virginia  22314 


Ooenanding  Officer 

Office  of  Naval  Research  Branch  Office 
$38  South  Clark  Street 
Chicago.  Illinot«  6060$ 


San  Francisco  Area  Office 
Office  of  Naval  Research 
$0  Fall  Street 

San  Francisco.  California  94102 


CoMkandlng  Officer 

Office  of  Naval  Research  Branch  Office 
1030  east  Green  Street 
Pasadena,  California  91X01 


CoMasnding  Officer 

Office  of  Havel  Rssearch  Branch  Office 

495  Stover  Street 

Boston.  Nassachusetts  02210 


New  York  Area  Office 
Office  of  Naval  Rssearch 
11$  M'oadMy  5th  Floor 
NSW  York,  New  York  10003 


ODDIUB  Mvlsory  Group  on  Electron  Devices 

201  varick  Street 

New  York.  New  York  10014 


Naval  Air  Davelopaent  Center 
Attn?  TWchnical  Library 
Johnsville 

NankinSCer.  Pennsylvania  1B974 


Naval  Naapons  Center 
Chine  Lake,  caiifronia  9359$ 
Attn:  Technical  Library  (I  copy) 
Qjrfe  6010  <2  copr^ 


Naval  Rasearch  Laboratory 

Underwater  Bound  Reference  Division 

TNetaleal  Library 

P.O.  Bos  8337 

Orlando.  Florida  32806 


Navy  Underwater  sound  Laboratory 
Technical  Library 
Port  TrtMiMjIl 

New  London.  Connecticut  06320 


CoHundant,  Marina  Corpa 
Sciantlfic  Adviaor  (Cado  AS) 
Haahington.  D.C.  20380 


Naval  Ordnance  Station 

Technical  Library 

Indian  Head,  Maryland  20640 


Naval  Postgraduata  School 
ttontaray,  California  93940 
Attn:  Technical  Library  fl  copy) 

Elect.  Bngin.  Depart.  (1  copy) 


Naval  Hlasile  Center 
Technical  Library  (Code  $632.2) 
Point  Nagu.  California  93010 


Naval  Electronics  Laboratory  Center 
San  Diego.  California 
Attn:  Technical  Library  (1  copy) 
Code  2300  (1  copy) 

2600  <1  copy) 

4800  (1  copy) 


Naval  Undersea  Center 

Technical  Library 

San  Olego.  California  92132 


Naval  Weapons  Laboratory 
Technical  Lll>rary 
OahXgren,  Virginia  22448 


Naval  Ship  Research  and  Developaent  Center 
Central  Library  (Codes  L42  and  L43) 
Nashington,  O.C.  20007 


Naval  Surface  Weapons  Center 
White  Oak  laboratory 
Silver  Spring,  ttoryland  20910 
Attn:  Technical  Library  (1  copy) 
Code  200  (1  copy) 

212  (I  copy) 


Deputy  Chief  of  Navel  (^rations 
(Developaent)' 

Technical  Analysis  and  Advisory  Group 
(Code  NOP-077D) 

Mashington.  D.C.  20350 


Oneunder 

Naval  Air  Systems  Coemand 
Mashington.  D.C. 

ATW:  Code  310  (I  copy) 
360  (1  c<N>y) 


Oowmander 

Naval  Electrontce  Syatama  Co— end 
Washington.  D.C.  20360 
Attn:  Code  304  (1  copy) 

310  (1  copy) 


Coa— nder 

Naval  Sea  SysteeM  Ooss Mind 
Naahingcon,  D.C.  20360 


Naval  Surface  Neepons  Center 

Attn:  Library 

Dehlgren.  Virginia  22448 


Air  Force  Office  of  Scientific  Reaearch 
Attn:  Electronic  and  Solid  State 
Sciences  Division 
Oepertasnt  of  the  Air  Force 
Nsshlngton.  O.C.  20333 


Air  Force  Nespon  Ukboratory 
TecAnient  li^ory 
Kirtlend  Air  Force  Base 
Albuguergue.  Mew  Momico  67117 


Air  Force  Avionics  Laboratory 
Air  Force  syatesie  Co—mnd 
Tachnlcal  Library 
Mright'Pattereon  Air  Force  Base 
Dayton,  Ohio  45433 


Air  Force  Cambridge  Reaearch  Laboratory 
L.C.  Hanecom  Field 
Technical  Library 
Cambridge,  Naaaachueette  02138 


Harry  DiasK>nd  Laboracories 
Technical  Library 

Connecticut  Avenue  at  Van  Ness,  M.H. 
Mashington.  D.C.  20438 


U.S.  Army  Reaearch  Office 
Box  CM.  Duke  Station 
Durham,  North  Carolina  27706 


Director 

U.S.  Army  Engineering  Research 
end  Development  Laboratories 
Port  Belvoir.  Virginia  22060 
Attn:  Technical  Dociaeents  Center 


Director  National  Bureau  of  standards 
Attn:  Technical  Library 
Mashington,  D.C.  20234 


Naval  Research  laboratory 
4555  Overlook  Avenue.  8.N. 
Mashington,  O.c.  2037$ 
Attn:  code  5300  (X  copy) 
7100  <1  copy) 
7900  (1  copy) 


Naval  Electronics  Laboratory  Center 
San  Diego.  California  92152 
Attn;  code  2100  (1  copy) 

2200  <1  copy) 


C.C.  Kliek 
Superintendent 
Hatariela  Sciences  Division 
Naval  Research  Laboratory 
4555  Overlook  Avenue.  S.M. 
Washington.  t>.C.  2037$ 


Navel  Research  Laboratory 
4555  Cberlook  Avenue.  S.M. 
Mashington.  D.C.  20375 
Attn:  Code  5220  (1  copy) 
$230  (1  copy) 
$250  (1  copy) 
$260  (1  copy) 
5270  (1  copy) 
5500  (1  ^py) 


Naval  Slsctronica  Laboratory  Canter 
San  Diego.  California  92152 
Attn;  Coda  2500  (1  copy) 

4000  (I  cofnr) 


Office  of  Naval  Reaearch  12  epa) 
800  N.  Quincy  street 
Arlington.  Virginia  22217 
Attn:  Coda  430  (2  copiaa) 


Nave)  Research  Laboratory 
4555  Overlook  Avenue,  S.M. 
Mashington.  D.C.  2037$ 
Attni  Code  5400 


Naval  Electronics  Latwratory  Canter 
San  Diego.  California  92152 
Attn:  Code  3000  (1  copy) 

5000  (I  copy) 

5600  (1  copy) 


Air  Force  office  of  Scientific  Reeaerch 
Natheamtlcal  and  information  Belancas 
Directorate 
1400  Nilaon  Blvd. 

Washington.  D.C.  20333 


